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The  Duroose  of  this  work  was  to  (i)  examine  separation  shock  wave  unsteadiness  in  different  turbulent  interactions  and  deter¬ 
mine  wLSJf  a  universi  model  describing  the  unsteadiness  could  be  developed,  and  (ii)  determine  whether  or  not  the  observed  un¬ 
steadiness  is  a  feature  of  turbulent  Bow  in  general,  or  is  specific  to  the  wind  tunnel  environment.  To  this  end.  wall  and  pitot  pressure 
fluctuation  measurements  were  made  in  interactions  generated  by  unswept  and  25  deg  swept  compression  ramp  models,  and  by  8  de„ 
and  30  de<J  sweet  blunt-fin  models  in  a  high  Reynolds  number.  Mach  5  turbulent  boundary  layer. 

It  is  clear  fhat  the  high-frequency,  jittery  motion  of  the  separation  shock  is  the  result  of  the  passage  through  the  wave  of  indi¬ 
vidual  larc^e-scale  turbulent  structures.  Thus,  this  component  of  the  unsteadiness  is  an  inherent  feature  of  afl  turbu^lent  flows  The 
primary  outstanding  question  concerns  the  cause  of  the  low-frequency  expansion/contraction  of  the  separated  flow  which  is  character¬ 
ized  bv  the  lar<je-scale.  long-duration  excursions  of  the  separation  shock  wave.  lo 

Preliminary  experimental  work  to  address  this  question  has  revealed  two  very  interesting,  complementary  results_  First,  there  s 
a  distinci  correction  between  large-scale  expansion  or  contraction  of  the  separated  flow  and  long  duration 

or  rises  in  pitot  pressure  in  the  incoming  turbulent  boundary  layer.  Second,  results  from  the  same  experiment  show  that  the  ensem- 
j  ble-averaged  pitot  pressure  at  a  fixed  location  in  the  incoming  undisturbed  boundary  layer  correlates  with  separation  shock  wave 

'^°^'The‘Je  results  suggest  a  model  in  which  the  boundary  layer  “thickens  and  thins"  as  it  flows  downstream.  However  at  this  siap  it 
must  be  emphasized  that  the  “thickening-thinning”  boundary-layer  model  is  an  inference  and  that  there  is  no  direct  time-resolv 
TviJence  thS  context  it  should  be  noted  that  some  evidence  has  been  found  of  Taylor-Gortler  vortices  in  the  incoming  flow  and 
it  is  possible  that  these  play  a  role  in  the  physics.  Confirmation  is  extremely  difficult  using  intrusive  instrurnentation  alone,  and.  in 
the  nCxt'  phase  of  the  work,  both  conventional  instrumentation  and  planar  laser  irnaging  techniques  will  be  used  to  explore  the 

incomino  flow  structure  corresponding  to  different  separated  bubble  states  and  dynamics. 
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Foreword 


The  nurnose  of  this  work  was  to  (i)  examine  separation  shock  wave  unsteadiness  in  different 
JiLtenuLSon*  and  deternUne  whether  a  universal  model  describing  the  unsteatoess 

could  be  developed,  and  (ii)  determine  whether  or  not  the  obsereed  unste^uess  is  a 

turbulent  flow  in  general,  or  is  specific  to  the  wind  tunnel  environment.  To  this  end,  wall  an 
pitot  pressure  fluctuation  measurements  were  made  in  interactions  generated  by  unswept  ^ 

S  deg  swem  compression  ramp  models,  and  by  8  deg  and  30  deg  swept  blunt-fm  inodels  m  a 
high  ReynoWs  number,  Mach  5  turbulent  boundary  layer.  Using  a  new  conditional  sampling 
algorithm,  the  separation  shock  unsteadiness  was  characterized  in  terms  of  its  and 

velocity  histories^  These  fluctuations  were  then  correlated  with  conditionally  extracted  static 
nressure  ratio  across  the  separation  shock,  and  with  wall  pressure  measurements  made  up- 
Leam  and  downstream  of  the  region  of  shock  unsteadiness.  Results  show  that  separation 
shock  unsteadiness  can  be  caused  by  two  different  physical  phenomena,  small-scale  or 
jittery  motion  is  caused  by  perturbations  in  the  ratio  of  static  quantities  across  the  shock  foot 
These  perturbations  induce  fluctuations  in  the  separation  shock  velocity  and  thus  changes  in 
lepTatSroc^  are  an  integral  result  of  the  velocity  fluctuations.  The  large-scale  or 
global  motion  of  die  separation  shock  is  caused  by  the  expansion  and  contraction  or  trembling 
S"  sepLion  bubble.  An  expansion  of  the  separation  bubble  or  exampi^  4  - 
nlaces  the  separation  shock  to  an  upstream  position,  whereas  a  contraction  of  the  bubble  result 
Fn  a  displacement  in  the  downstream  direction.  In  addition  to  experimental  results,  ^omp  - 
tions  of  an  inviscid,  unsteady,  quasi  one-dimensional  flow  have  also  been  performed. 
Although  this  does  not  model  the  experimental  conditions,  the  analysis  provides  ^ompmationa 
support  for  the  explanation  offered  for  the  small-scale  jittery  motion  of  the  separation  shock 

As  described  above,  it  is  clear  that  the  high-frequency,  jittery  motion  of  the  separation 
shock  is  the  result  of  the  passage  through  the  wave  of  individual  large-scale  turbulent  stru  - 
tures  Thus,  it  is  reasonable  to  conclude  that  this  component  of  the  unsteadiness  jnheren 
fpatiire  of  all  turbulent  flows.  The  primary  outstanding  question  concerns  the  cause  of  the  low 

of  the^rated  flow  which  is  characterized  by  the  large- 

scale,  long-duration  excursions  of  the  separation  shock  wave.  •  , 

Preliminary  experimental  work  to  address  this  question  has  revealed  two  very  interest- 
ma  rnmnlernentarv  results.  First,  there  is  a  distinct  correlation  between  large-scale  expansion 
or^contraction  of  the  separated  flow  and  long  duration  (i.e.,  low-frequency)  falls  or  ^ses  m 
nitot  nressure  in  the  incoming  turbulent  boundary  layer.  The  time  scale  of  these  gradual  ns 
Ld  falls  in  pitot  pressure  is  around  1  ms,  almost  two  orders  of  magnitude  higher  than  *e  large 
eddv  time  scale  5 1  lUoc  Thus,  these  rises  and  falls  in  pitot  pressure  are  not  generate  Y 
vM^rbufut  st^^^^^^  but  are  caused  by  low-frequency  “global”  motion  of  the  bound^ 
laver  Second  results  from  the  same  experiment  show  that  the  ensemble-averaged  pi  p 
s^a.  afeSSta  [n  fte  incoming  undisturbed  boundary  layer  correlates  w.th  separatton 

shock  a  model  in  which  the  boundary  layer  “Ihickens  and 

thins”  as  it  flows  downstream.  In  such  a  simple  model,  a  higher  pitot  pressure  at  a  fixed  heig 
w  JildTrtlply  a  X^ity  profile  than  the  mean  profile,  one  more  resistant  to  separation. 
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which  would  generate  a  smaller  separation  bubble  than  would  the  mean  profile.  Conversely,  a 
lower  pitot  pressure  at  a  fixed  height  would  imply  a  slightly  retarded  profile  which  would  gen¬ 
erate  a  larger  separation  bubble  than  would  the  mean  profile.  However,  at  this  stage  it  must  be 
emphasized  that  the  “thickening-thinning”  boundary-layer  model  is  an  inference  and  that  there 
is  no  direct  time-resolved  evidence.  In  this  context  it  should  be  noted  that  some  evidence  has 
been  found  of  Taylor-Gortler  vortices  in  the  incoming  flow,  and  it  is  possible  that  these  play  a 
role  in  the  physics.  Confirmation  is  extremely  difficult  using  intrusive  instrumentation  alone, 
and,  in  the  next  phase  of  the  work,  both  conventional  instrumentation  and  planar  laser  imaging 
techniques  will  be  used  to  explore  the  incoming  flow  strucmre  corresponding  to  different  sepa¬ 
rated  bubble  states  and  dynamics. 
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Nomenclature 


Cf  skin  friction  coefficient 

D  cylinder  or  blunt-fin  leading  edge  diameter 
G  power  spectral  density  coefficient 
H  shape  factor 

M  Mach  number 

P  pressure 

R  cross-correlation  coefficient 

Rc  computed  pressure  ratio  across  normal  shock 

Rs  measured  pressure  ratio  across  separation  shock 

Re  unit  Reynolds  number 

T  threshold 

U  velocity 

Vc  computed  normal  shock  velocity 

Vs  measured  separation  shock  velocity 

Xc  computed  normal  shock  position 

Xs  measured  separation  shock  position 

f  fall-time  or  frequency 

h  channel  height 

p  pressure 

r  rise-time 

t  time 

X  streamwise  distance 

z  spanwise  distance 

AY  ordinate  offset 

n  boundary  layer  wake  strength  parameter 

5o  boundary  layer  thickness 

5*  boundary  layer  displacement  thickness 

X  sweepback  angle 

Y  separation  shock  intermittency  factor 

0  boundary  layer  momentum  deficit  thickness 

p  density 


m 


a  standard  deviation 

X  time-delay 
C  transducer  spacing 

(  )  mean  value 

Subscripts 

EA  ensemble-averaged 
d  downstream 

i  channel  inlet 

0  channel  outlet  or  incoming  conditions 

t  stagnation  quantities 
u  upstream 

w  wall  conditions 

oo  freestream  conditions 

Abbreviations 

Ch  channel 

FIR  Finite  Impulse  Response 

LDV  Laser  Doppler  Velocimeter 

S  separation 

SBF-30  30  deg  swept  blunt-fin 

SBF-8  8  deg  swept  blunt-fm 

SCR-25  25  deg  swept  compression  ramp 

UCR  unswept  compression  ramp 

UI  upstream  influence 

VITA  Variable  Interval  Time  Averaging 


IV 


1.  Introduction 


Shock  wave /turbulent  boundary  layer  interactions  are  important  in  a  wide  variety  of 
internal  and  external  aerothermodynamics  problems  of  interest  to  the  U.S.  Army.  They  include 
flows  over  transonic  rotor  blades,  around  high  speed  missile  control  surfaces  and  fins,  and 
upstream  of  a  missile  base  where  unsteady  separation  is  induced  by  the  plume.  Plume-induced 
separation  effects  can  play  a  critical  role  in  determining  missile  base  drag  and  plume  signature. 
Consequently,  there  has  been  a  very  large  number  of  studies  of  these  interactions  over  the  past 
four  decades.  Although  high  speed  cinematography  from  the  1950's  (Bogdonoff  1955, 
Chapman  et  al  1958)  showed  that  the  shock-induced  turbulent  separation  was  unsteady,  the 
lack  of  adequate  instrumentation  resulted  in  most  studies  being  of  mean  flow  properties  only. 
Kistler  (1964)  was  probably  the  first  to  make  detailed  wall  pressure  measurements  under  sepa¬ 
rated  supersonic  turbulent  boundary  layer,  and  the  first  to  discuss  the  intermittent  character  of 
the  wall  pressure  signal  caused  by  the  translating  separation  shock  wave.  Since  Kistler  there 
has  been  an  increasing  number  of  studies  focusing  on  interaction  unsteadiness.  The  qualitative 
character  and  quantitative  details  of  the  unsteadiness  naturally  vary  from  one  flow  type  to  an¬ 
other,  but  typically  the  unsteadiness  manifests  itself  as  a  large-scale,  low-frequency  pulsation 
of  the  separated  flow  and  flapping  of  the  outgoing  boundary  layer.  A  recent  review  of  much  of 
this  work  has  been  compiled  by  Dolling  (1993). 

From  an  engineering  perspective  the  unsteadiness  is  important  from  two  points  of 
view.  First,  it  results  in  large  amplitude  fluctuating  loads  (as  high  as  185  dB)  which  can  sub¬ 
stantially  shorten  the  fatigue  life  of  vehicle  components.  Second,  it  compounds  the  difficulties 
of  drawing  conclusions  from  comparisons  of  experimental  data  with  computations.  The  for¬ 
mer  are  largely  time-averaged  data  from  an  unsteady  flowfield  whereas  the  latter  represent  con¬ 
verged  or  steady  state  flowfields.  The  difficulties  of  drawing  meaningful  conclusions  from 
such  comparisons  are  discussed  by  Dolling  (1992)  and  Marshall  &  Dolling  (1992).  Clearly,  an 
understanding  of  the  physical  causes  of  the  unsteadiness  is  needed,  not  only  to  understand  the 
flowfield  physics  but  to  assist  in  the  development  of  appropriate  modeling  strategies.  At  the 
present  time,  understanding  of  the  physical  mechanisms  is  incomplete,  although  as  outlined 
below  there  has  been  progress  over  the  past  few  years. 

Andreopoulos  &  Muck  (1987)  made  wall  pressure  measurements  in  Mach  3  compres¬ 
sion  ramp  interactions,  and  reported  that  the  mean  separation  shock  wave  period  was  indepen¬ 
dent  of  position  in  the  intermittent  region  and  independent  of  the  ramp  angle.  The  intermittent 
region  is  the  region  within  which  the  separation  shock  translates.  The  zero-crossing  fre¬ 
quency,  which  is  simply  the  number  of  unidirectional  shock  crossings  over  a  given  transducer 
in  the  intermittent  region,  was  found  to  be  0.13  Uoo/5o,  where  Uoo  is  the  freestream  velocity 
and  5o  is  the  incoming  boundary  layer  thickness.  Since  this  frequency  was  the  same  order  as 
the  estimated  bursting  frequency  in  the  incoming  turbulent  boundary  layer  and  since  the  mea¬ 
sured  shock  speeds  were  the  same  order  as  velocity  fluctuations  in  the  flowfield,  Andreopoulos 
&  Muck  argued  that  “the  incoming  boundary  layer  is  largely  responsible  for  the  shock 
motion.”  If  the  shock  frequency  is  controlled  by  turbulent  bursts  and/or  velocity  fluctuations, 
then  this  reasoning  would  suggest  that  in  a  given  incoming  boundary  layer  the  shock  zero¬ 
crossing  frequency  and  the  streamwise  extent  of  its  motion  would  be  fixed.  This  certainly  is 
not  the  case.  For  cylinder-induced  interactions  in  a  fixed  incoming  boundary  layer,  the  zero¬ 
crossing  frequency  decreases  and  the  streamwise  extent  of  the  shock  motion  increases  as  the 
cylinder  diameter  D  increases  (Dolling  &  Smith  1989).  In  fact,  the  intermittent  region  length  is 
essentially  0.8  D,  independent  of  boundary  layer  thickness.  Furthermore,  since  the  one- 
threshold  technique  used  by  Andreopoulos  &  Muck  has  been  shown  to  overestimate  shock  fre¬ 
quencies  significantly  because  of  its  inability  to  discriminate  between  shock-induced  and  turbu¬ 
lent  pressure  fluctuations  (Dolling  &  Brusniak  1989),  it  is  probable  that  changes  in  shock 
frequency  with  ramp  angle  were  not  detected. 
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In  a  later  study,  Tran  (1986)  made  simultaneous  wall  pressure  measurements  under  the 
incoming  boundary  layer  and  near  the  upstream  influence  line  in  the  interaction  generated  by  a 
20  deg  compression  ramp  model.  The  experiments  were  made  in  the  same  test  facility  used  by 
Andreopoulos  &  Muck  (1987).  Using  the  Variable  Interval  Time  Averaging  (VITA)  technique, 
Tran  found  little  correlation  between  “events”  detected  on  the  upstream  channel  and  the  shock- 
induced  pressure  pulses  on  the  downstream  channel.  He  concluded  that  the  pressure  pulses  in 
the  intermittent  region  were  independent  of  the  large-scale  turbulent  structures  in  the  upstream 
boundary  layer  which  are  convected  into  the  interaction.  Tran's  conclusions,  however,  must 
also  be  treated  with  some  caution.  He  used  only  a  single  transducer  in  the  intermittent  region, 
close  to  its  upstream  edge.  It  is  probable  that  the  separation  shock  moved  upstream  and 
downstream  in  response  to  the  turbulent  structures,  but  remained  downstream  of  the  measuring 
transducer.  With  hindsight,  it  would  be  better  to  use  the  shock-induced  pressure  pulse  on 
downstream  transducer  as  the  trigger  and  ensemble  average  on  the  upstream  channel  in  the 
undisturbed  boundary  layer. 

This  approach  was  adopted  by  Erengil  &  Dolling  (1991b)  and  later  refined  by  McClure 
(1992).  Simultaneous  measurements  were  made  of  wall  pressure  fluctuations  under  the  incom¬ 
ing  boundary  layer  and  the  intermittent  region  of  a  Mach  5  compression  ramp  interaction. 
Conditional  sampling  algorithms  and  a  variable  window  ensemble  averaging  technique  were 
used  to  investigate  the  correlation  between  specific  separation  shock  motions  and  pressure 
fluctuations  in  the  incoming  turbulent  boundary  layer.  Although  the  correlation  was  clearly 
evident  for  turn-around  motions  of  the  separation  shock,  the  interpretation  of  the  results  corre¬ 
sponding  to  upstream  and  downstream  sweeps  was  not  clear  due  to  an  insufficient  number  of 
ensembles.  Nevertheless,  to  the  best  of  the  authors’  knowledge,  these  results  were  probably 
the  first  to  show  a  direct  correlation  between  pressure  fluctuations  in  the  incoming  boundary 
layer  and  the  separation  shock  motion. 

The  interpretation  of  the  results  corresponding  to  upstream  and  downstream  motions  of 
the  separation  shock  has  since  been  clarified  by  McClure  (1992).  The  experimental  arrange¬ 
ment  was  the  same  as  that  of  Erengil  &  Dolling  (1991b).  McClure  combined  sweep  motions 
from  different  positions  in  the  intermittent  region,  creating  a  large  number  of  ensembles  for  a 
given  sweep  motion,  say  an  upstream  sweep  over  three  transducers.  These  results  showed  that 
the  correlation  seen  in  the  turn-around  motions  was  not  unique  to  these  motions,  and  that  simi¬ 
lar  correlations  also  exist  for  upstream  and  downstream  sweep  motions  of  the  shock.  Most 
importantly,  the  characteristic  shapes  of  the  ensemble  signatures  on  the  upstream  channels  are 
different  for  the  different  shock  motions.  Specifically,  the  pressure  signatures  for  the  upstream 
sweeps  have  a  rise-fall-rise  sequence,  whereas  those  for  the  downstream  sweeps  have  a  fall- 
rise-fall  sequence.  Also,  analysis  of  the  timing  of  events  showed  that  these  signatures  are  con¬ 
vected  into  the  interaction  and  are  coincident  with  the  separation  shock  foot  at  the  time  of  the 
crossing  over  a  given  transducer.  The  structures  causing  these  signatures  convect  downstream 
pass  through  the  shock  and  continue  downstream  through  the  interaction  at  convection  speeds 
of  about  0.75  Uoo.  Their  temporal  extent  is  about  75  to  100  ps,  which  corresponds  to  a  spa¬ 
tial  extent  of  about  3.2  to  4.3  5o.  In  an  extension  of  McClure’s  work,  Gramann  &  Dolling 
(1992)  detected  the  signatures  as  far  upstream  as  20  5o  from  the  interaction. 

Gramann  (1989)  studied  the  dynamics  of  the  separation  bubble  in  a  Mach  5  compres¬ 
sion  ramp  interaction.  Using  conditional  sampling  algorithms  as  well  as  standard  analysis 
techniques  applied  to  wall  pressure  fluctuation  measurements,  he  specifically  examined  the 
relative  motion  of  the  separation  and  reattachment  lines  in  this  interaction.  His  results  showed 
that  the  separation  bubble  undergoes  an  expansion  and  contraction  type  motion  in  which  the 
separation  and  reattachment  lines  move  away  from  and  towards  each  other.  Similar  observa¬ 
tions  were  also  made  earlier  by  Kussoy  et  al.  (1987)  in  a  Mach  3  cylinder-flare  interaction.  In 
the  latter  study,  high  speed  shadow  movies  and  LDV  measurements  were  made  simultaneously 
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with  wall  pressure  measurements  under  the  unsteady  separation  shock.  Results  on  the  plane  of 
symmetry  showed  that  the  separation  bubble  “expanded  and  contracted  like  a  balloon.” 

In  another  study,  Selig  &  Smits  (1991)  applied  periodic  blowing  to  the  boundary  layer 
at  various  locations  within  the  separated  region  generated  by  a  24  deg  compression  ramp  in  a 
Mach  3  flow.  Blowing  was  applied  through  a  slot  across  the  full  span  at  a  rate  of  either  2.5  or 
9  percent  of  the  freestream  mass  flux  at  frequencies  up  to  5  kHz.  Wall  pressure  and  mass  flux 
fluctuation  measurements  were  made  and  stroboscopic  schlieren  videography  was  used  to  ob¬ 
serve  the  unsteady  flowfield.  As  the  slot  was  moved  closer  to  the  separation  shock  foot,  the 
shock  motion  could  be  locked  in  to  the  blowing  frequency  and  completely  controlled  by  the 
periodic  blowing.  Although  the  upstream  influence  had  increased,  the  rms  levels  in  the  vicinity 
of  the  separation  shock  did  not  change  appreciably.  Based  on  cross-correlation  results  and 
their  success  in  the  control  of  the  unsteadiness,  Selig  and  Smits  concluded  that  “it  seems  prob¬ 
able  that  the  unsteadiness  of  the  shock  for  the  24  deg  compression  comer  is  partially  driven  by 
large-scale  fluctuations  in  the  separated  region.” 

Although  some  of  these  earlier  studies  show  a  direct  correlation  between  separation 
shock  motion  and  pressure  fluctuations  in  the  incoming  turbulent  boundary  layer,  they  fail  to 
provide  a  physical  model  by  which  these  turbulent  stmctures  are  related  to  the  separation  shock 
unsteadiness.  By  the  same  token,  although  those  studies  investigating  a  possible  correlation 
between  separation  shock  motion  and  dynamics  of  the  separation  bubble  have  demonstrated  the 
characteristic  motion  of  the  latter,  they  too  have  not  been  successful  in  providing  a  physical 
model  which  explains  the  unsteadiness.  This  is  a  particularly  difficult  task  since  separation 
shock  motion  and  dynamics  of  the  separation  bubble  as  well  as  fluctuations  in  the  separated 
flow  are  highly  coupled. 

The  purpose  of  the  study  reported  in  this  paper,  is  to  examine  separation  shock 
unsteadiness  in  different  interactions  and  first  determine  whether  a  universal  model  describing 
and  explaining  the  unsteadiness  could  be  developed.  The  second  objective  was  to  try  and 
determine  if  the  results  are  an  inherent  feature  of  all  turbulent  flows  or  whether  they  are  a  fea¬ 
ture  seen  solely  in  wind-tunnel-generated  boundary  layers.  To  do  this,  existing  measurements 
from  different  interactions  have  been  used  and  also  new  experiments  performed  whenever  nec¬ 
essary.  Measurements  from  unswept  and  25  deg  swept  compression  ramp  flows,  and  from  8 
deg  and  30  deg  swept  blunt-fin  interactions  have  been  analyzed.  Using  simultaneously  sam¬ 
pled  multi-channel  wall  pressure  fluctuations,  separation  shock  foot  position  history,  Xs(t)  and 
velocity  history,  Vs(t)  have  been  determined.  Once  they  were  determined,  these  histories  were 
correlated  with  conditionally  extracted  static  pressure  ratio,  Rs(t),  across  the  separation  shock, 
and  with  wall  pressure  measurements  made  upstream  and  downstream  of  the  interrnittent 
region.  In  this  paper,  the  discussion  of  these  and  other  results  are  followed  by  the  description 
of  a  physical  model  describing  the  separation  shock  unsteadiness  in  shock  wave  /  turbulent 
boundary  layer  interactions. 

In  addition  to  the  experimental  results,  computations  of  an  unsteady,  inviscid  and  quap 
one-dimensional  flow  are  also  presented.  It  should  be  emphasized  that  the  purpose  of  this 
work  is  not  to  model  the  unsteadiness  in  these  interactions.  Instead,  the  CFD  calculations  were 
used  as  a  tool  to  further  demonstrate  one  part  of  the  physical  mechanism  of  unsteadiness 
obtained  from  experimental  results. 

2.  Experimental  Program 

2.1.  Wind  tunnel  and  flow  conditions 

The  experiments  were  conducted  in  the  high  Reynolds  number,  Mach  5  blowdown 
wind  tunnel  at  Balcones  Research  Center  of  the  University  of  Texas  at  Austin.  The  tunnel  has 
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a  rectangular  test  section  which  is  15.24  cm  (6  in.)  wide,  17.78  cm  (7  in.)  high,  and  30.48  cm 
(12  in.)  long.  Up  to  3.964  m^  (140  fP)  of  compressed  air  stored  at  about  17.24  MPa  (2500 
psia)  provides  run  times  of  up  to  1  min.  for  the  current  stagnation  conditions.  The  incoming 
air  is  heated  by  two  420  kW  banks  of  nichrome  wire  resistive  heaters  located  upstream  of  the 
stagnation  chamber.  These  heaters  can  provide  stagnation  temperatures  of  up  to  422  K  (759 

°R). 

All  tests  were  conducted  at  a  stagnation  pressure  of  2.28  MPa  (330  psia)  and  a  stagna¬ 
tion  temperature  of  359  K  (647  oR),  which  generate  a  freestream  unit  Reynolds  number  of 
49  9  X 10^  m‘^  (15.2  x  10^  ff^).  The  freestream  Mach  number  and  velocity  in  the  test  sec¬ 
tion  were  4.95  and  775  m/s  (2543  ft/s),  respectively.  The  incoming  turbulent  boundary  layer 
underwent  natural  transition  and  developed  under  approximately  adiabatic  wall  temperature 
conditions  Boundary  layer  parameters  were  determined  from  pitot  surveys  made  with  a  probe 
having  a  tip  orifice  0.020  cm  (0.008  in.)  high  and  0.038  cm  (0.015  in.)  wide.  Total  tempera¬ 
ture  surveys  were  also  made  using  thermocouple  probes.  In  the  data  reduction,  the  static  pres¬ 
sure  through  the  boundary  layer  was  assumed  constant  and  equal  to  the  freestream  value.  A 
good  fit  to  law-of-the-wall  /  law-of-the-wake  was  obtained  for  the  mean  velocity  profile. 
Boundary  layer  parameters  measured  at  96.47  cm  (37.98  in.)  downstream  of  the  throat  on  the 
tunnel  centerline,  approximately  3  5o  upstream  of  the  models,  are  listed  in  table  1  (McClure 

1992). 
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Cf 
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Table  1 .  Boundary  layer  parameters. 


2.2.  Models  and  transducer  plug 

Several  models  have  been  used  to  generate  different  interactions.  These  include 
unswept  and  25  deg  swept  compression  ramp  models,  and  8  deg  and  30  deg  swept  blunt-tin 

The  unswept  28  deg  compression  ramp  model  was  made  of  aluminum  and  was  12.07 
cm  (4  75  in )  wide  and  3.8  cm  (1.5  in.)  high.  Two  aluminum  aerodynamic  fences  with  sha^- 
ened  leading  edges  were  attached  to  the  sides  of  the  ramp  to  prevent  spillage  and  isolate  the 
interaction  from  the  sidewall  boundary  layers.  The  25  deg  swept 

was  made  of  brass  and  was  10.8  cm  (4.25  in.)  wide  and  5.1  cm  (2  in.)  Ingh.  It  had  a  constant 
streamwise  angle  of  28  deg.  For  the  swept  ramp  only  one  stainless  steel  fence  was  attached  to 
the  apex  side  of  the  model  and  the  fence  leading  edge  extended  at  least  6.4  cm  (2.5  in.)  up- 
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stream  of  the  apex.  Both  blunt-fm  models  were  made  of  aluminum.  They  had  a  leading  edge 
diameter  of  1.9  cm  (0.75  in.)  and  were  10.16  cm  (4  in.)  high.  Ail  models  were  mounted  on 
the  tunnel  floor  and  could  be  moved  in  the  streamwise  direction.  Sketches  of  these  models  and 
their  appropriate  coordinate  systems  are  shown  in  Fig.  1. 

The  brass  instrumentation  plug,  which  was  8.57  cm  (3.375  in.)  in  diameter,  was  cen¬ 
tered  97.00  cm  (38.19  in.)  downstream  of  the  nozzle  throat  and  installed  in  the  test  section 
floor.  In  this  plug,  the  transducers  could  be  mounted  flush  along  four  different  rows  which 
could  be  aligned  at  any  angle  to  the  freestream.  The  first  row  was  along  the  plug  centerline  and 
it  had  26  transducer  ports.  The  second  (8  ports),  third  (13  ports),  and  fourth  (12  ports)  were 
parallel  to  the  center  row  and  were  spaced  0.29  cm  (0.1 15  in.),  1.27  cm  (0.5  in.),  and  2.54  cm 
(1  in.)  from  the  latter,  respectively.  In  any  given  row,  the  center-to-center  spacing  between 
two  adjacent  ports,  C.  was  0.29  cm  (0.115  in.).  Customized  dummy  brass  plugs  were 
mounted  flush  in  the  unused  transducer  ports. 

2.3.  Instrumentation  and  data  acquisition  system 

Up  to  eight  miniature,  Kulite,  high-frequency  response  pressure  transducers  were 
used.  These  were  either  XCW-062-15A,  0-103  kPa  (0-15  psia)  or  XCQ-062-50A,  0-345  kPa 
(0-50  psia)  absolute  pressure  transducers.  Transducers  were  installed  either  flush  with  the 
floor  to  measure  fluctuating  wall  pressure  or  inside  of  small  probes  used  to  measure  fluctuating 
pitot  pressure.  They  all  had  outer  diameters  of  0.159  cm  (0.0625  in.)  and  0.071  cm  (0.028 
in.)  diameter  pressure  sensitive  diaphragms.  As  quoted  by  the  manufacturer,  the  natural  fre¬ 
quency  response  of  the  diaphragm  was  250  kHz  for  0-103  kPa  models  and  500  kHz  for  0-345 
kPa  models.  However,  perforated  screens  protecting  the  transducers  from  dust  particles  limit 
the  effective  frequency  response  of  both  models  to  about  50  kHz.  All  transducers  were  cali¬ 
brated  statically  using  a  Heise  digital  pressure  gauge  (Model  710A)  accurate  to  within  6.9  Pa 
(0.001  psia). 

Output  from  pressure  transducers  was  first  amplified  by  either  Dynamics  (Model 
7525),  Measurements  Group  (Model  2311),  or  PARC  (Model  113)  amplifiers,  and  then  low- 
pass  filtered  using  Ithaco  (Model  4113  or  Model  4213)  analog  filters.  Amplifier  gains  were 
adjusted  for  maximum  resolution  and  varied  depending  on  the  transducer  model  and  its  posi¬ 
tion  in  the  interaction.  Filter  cut-off  frequencies  were  set  at  50  kHz  unless  the  sampling 
frequency  required  a  lower  setting. 

Filtered  signals  were  then  digitized  by  two  LeCroy  Model  6810  Waveform  Recorders 
equipped  with  12-bit  analog-to-digital  (A/D)  converters.  Each  of  these  recorders  had  two 
Megawords  of  memory,  and  could  sample  1  channel  of  data  at  rates  up  to  5  MHz  per  channel, 
or  4  channels  of  data  simultaneously,  at  rates  up  to  1  MHz  per  channel.  Both  A/D  converters 
were  triggered  using  the  same  clock,  thus  either  2  channels  could  be  sampled  simultaneously  at 
rates  up  to  5  MHz  per  channel,  or  8  channels  at  rates  up  to  1  MHz  per  channel.  A  maximum  of 
4096  records  (1  record  =  1024  datapoints)  of  data  could  be  obtained.  Further,  each  channel 
could  be  adjusted  independently  allowing  the  user  to  optimize  each  channel  for  different  input 
signal  ranges. 

Data  were  then  downloaded  to  an  HP  9000  Model  380  workstation  and  stored  on  mag¬ 
netic  tapes  for  later  analysis.  All  subsequent  analyses  were  performed  either  on  the  primary 
computer  (Model  380)  or  on  one  of  the  two  HP  9000  Model  340  workstations  connected  to  the 
primary  computer. 


5 


2.4.  Test  program 

The  main  purpose  of  this  work  was  to  examine  the  separation  shock  unsteadiness  ob¬ 
served  in  different  interactions,  to  determine  whether  a  universal  model  describing  the  causes 
of  unsteadiness  could  be  developed,  and  to  understand  the  source  of  the  unsteadiness.  To  this 
end,  the  authors  made  use  of  existing  measurements  performed  by  themselves  and  their  col¬ 
leagues  in  the  same  facility.  Experiments  in  the  unswept  and  swept  compression  ramp  interac¬ 
tions  were  performed  by  Erengil  &  Dolling  (1993).  Those  in  swept  blunt-fm  interactions  were 
done  by  Kleifges  &  Dolling  (1993).  In  addition  to  these,  some  new  experiments  have  been 
performed  by  the  authors  in  interactions  generated  by  the  unswept  compression  ramp  and  the  8 
deg  swept  blunt-fm  models.  In  all  these  experiments,  at  least  six  transducers  were  placed  up¬ 
stream  of  the  separation  line  to  measure  the  separation  shock  dynamics.  The  location  of  the 
separation  line  (i.e.,  line  of  coalescence  of  surface  streaklines)  is  usually  determined  from  sur¬ 
face  tracer  techniques.  In  this  case,  kerosene-diesel-lampblack  mixture  is  painted  on  the  tunnel 
floor  prior  to  a  flow  visualization  mn,  and  the  pattern  is  lifted  off  the  surface  on  large  sheets  of 
transparent  tape.  The  remaining  transducers  were  used  to  measure  wall  pressure  fluctuations 
under  the  incoming  turbulent  boundary  layer  and/or  under  the  separated  shear  layer. 

In  all  the  interactions  studied,  the  unsteadiness  is  characterized  by  the  separation  shock 
position,  Xs(t),  and  velocity,  Vs(t),  histories  which  are  derived  from  simultaneous  ineasure- 
ments  of  wall  pressure  flucmations  in  the  intermittent  region.  Consequently,  to  determine  Xs(t) 
and  Vs(t)  accurately,  one  must  use  as  many  transducers  as  possible  within  the  intermittent  re¬ 
gion  of  a  given  interaction.  This  will  be  clearer  as  the  derivation  of  these  quantities  is  explained 
in  the  following  section. 

3.  Analysis  Techniques 

Standard  time-series  analysis  techniques  as  well  as  conditional  sampling  algorithms 
have  been  used  to  analyze  the  experimental  data.  In  the  conditional  sampling  techniques  em¬ 
ployed  here  the  separation  shock  position,  Xs(t),  and  velocity,  Vs(t),  histories  are  calculated. 
Before  discussing  the  derivation  of  these  quantities  in  some  detail,  it  will  be  helpful  to  first  de¬ 
scribe  the  box-car  conversion  technique  briefly  since  Xs(t)  and  Vs(t)  are  deduced  from  nested 
box-car  sequences.  A  detailed  description  of  the  two-threshold  box-car  conversion  technique 
has  been  published  by  Dolling  &  Bmsniak  (1989). 

3.1.  Box-car  conversion  technique 

The  two-threshold  box-car  conversion  technique  is  used  to  extract  shock  passages  from 
wall  pressure  fluctuations  measured  under  the  unsteady  separation  shock.  In  this  method,  the 
thresholds  are  set  at  Ti  =  Pwo  +  3(7t’wo  and  Tz  =  Pwo  +  60lPwo,  where  Pwo  is  the  undisturbed 
turbulent  boundary  layer  mean  pressure  level  and  OTPwo  is  the  rms  of  the  turbulent  boundary 
layer  pressure  fluctuations.  These  threshold  settings  were  recommended  following  a  sensitiv¬ 
ity  analysis  of  the  effects  of  varying  Ti  and  Tz  (Dolling  &  Brusniak  1989).  A  siinilar 
sensitivity  analysis  showed  that  the  same  threshold  settings  were  appropriate  for  the  analysis  of 
data  from  swept  flows  (Erengil  &  Dolling  1993).  In  the  two-threshold  method,  the  box-car  is 
generated  when  Pw(t)>Tz  and  is  terminated  when  Pw(t)<Ti.  Thus  when  the  box-car 
assumes  a  value  of  1,  this  indicates  a  shock  passage  in  the  upstream  direction  and  is  called  the 
rise-time.  Similarly,  when  it  assumes  a  value  of  0,  this  indicates  a  shock  passage  in  the 
downstream  direction  and  is  called  the  fall-time.  Once  a  shock  passage  is  detected  using  Ti 
and  Tz  and  an  initial  fall-time  assigned,  a  counter  is  then  used  to  march  backward  in  time  until 
the  instantaneous  value  of  Pw(t), which  is  just  below  Tz,  is  found.  The  fall-time  is  then  reas¬ 
signed  to  this  new  value.  The  requirement  that  the  rise-  and  fall-times  be  determined  by  the 
same  threshold  level  (i.e.,  Tz)  ensures  that  the  separation  shock  foot  is  at  the  same  position  at 
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both  rise-  and  fall-times.  The  relative  position  of  the  transducer  and  its  rise-  and  fall-times  are 
the  only  pertinent  information  for  the  derivation  of  the  separation  shock  position  history. 
Typical  wall  pressure  signals  from  the  unswept  compression  ramp  interaction  and  correspond¬ 
ing  box-car  sequences  are  shown  in  Figs.  2a  through  2f. 


3.2.  Derivation  of  separation  shock  position  and  velocity  histories 


Separation  shock  position  history  is  deduced  from  nested  box-car  sequences  which  are 
generated  from  simultaneous  measurements  of  multi-channel  wall  pressure  fluctuations  made  in 
the  intermittent  region.  The  schematic  at  top  of  Fig.  2  shows  the  position  of  eight  transducers, 
spaced  1^  apart,  upstream  of  the  separation  line  in  the  interaction  generated  by  the  unswept 
compression  comer.  Given  eight  transducers,  there  are  nine  discrete  bins  in  which  the  separa¬ 
tion  shock  foot  can  be  bracketed.  In  this  case,  a  bin  is  defined  as  the  spatial  region  from  the 
downstream  end  of  one  transducer  to  the  downstream  end  of  the  adjacent  transducer  as  indi¬ 
cated  in  the  schematic.  Simultaneous  wall  pressure  fluctuations  and  corresponding  box-car 
sequences  measured  at  stations  indicated  by  filled  markers  are  plotted  in  Figs.  2a  through  2f  for 
channels  1  through  6,  respectively.  As  was  discussed  earlier,  the  rise-  and  fall-times  of  a  given 
box-car  sequence  indicate  shock  passages  in  the  upstream  and  downstream  directions.  Some 
of  these  times  are  indicated  in  Figs.  2c  through  2e.  The  box-car  sequences  shown  in  Figs.  2a 
through  2f  are  plotted  again  in  the  nested  form  in  Fig.  2g.  The  rise-  and  fall-times  in  the  nested 
box-car  sequence  and  the  channel  on  which  they  occur  define  the  time-  and  bin-coordinates  of 
the  separation  shock  foot.  At  times  t  =  r5  and  t  =  r4,  for  example,  the  separation  shock  enters 
bins  5  and  4  from  bins  6  and  5,  respectively.  This  information  provides  two  points  in  the 
separation  shock  foot  history,  Xs(t)  —  indicated  by  rs  and  r4  in  Fig.  2h.  Due  to  the  spatial 
resolution  restrictions,  the  separation  shock  motion  cannot  be  determined  for  rs  <  t  <  r4  .  All 
that  is  known  is  that  it  stays  in  bin  5.  It  is  assumed  that  the  motion  of  the  separation  shock 
during  this  time  is  unidirectional  and  at  constant  speed.  This  assumption  allows  for  the  linear 
interpolation  of  the  two  time/bin  coordinates,  rs  and  r4 ,  as  shown  in  Fig.  2h.  Similarly,  an¬ 
other  interpolation  is  performed  from  t  =  r4  to  t  =  r3  (i.e.,  the  time  when  separation  shock 
enters  bin  3  from  bin  4).  At  time  t  =  f3  the  separation  shock  enters  bin  4  from  bin  3.  In  this 
case,  a  different  linear  interpolation  is  required  since  the  separation  shock  changes  its  direction 
of  motion  within  this  time  interval  (i.e.,  r3  <  t  <  f3).  It  is  further  assumed  that  the  separation 
shock  moves  to  the  middle  of  the  bin,  changes  its  direction  of  motion  and  moves  back  at  the 
same  speed.  This  assumption  results  in  the  interpolation  depicted  in  Fig.  2h.  Those  motions 
opposite  of  the  ones  discussed  above  (i.e.,  a  sequence  of  fall-times  and  a  fall-rise  pair)  ^e 
treated  in  a  similar  manner.  Clearly,  the  linear  interpolation  yields  a  piecewise  smooth  function 
for  Xs(t)  as  shown  in  Fig.  2h.  It  should  be  noted  that  the  assumptions  invoked  in  this  deriva¬ 
tion  are  not  new  and  have  always  been  implicit  in  conditional  analysis  techniques  employed 


previously.  ,  .  ,  ,  ... 

Once  the  separation  shock  position  history,  Xs(t),  is  obtained,  the  separation  shock 
velocity  history,  Vs(t),  can  easily  be  calculated  by  taking  the  derivative  of  Xs(t).  The  authors 
are  aware  of  the  problems  involved  in  differentiating  the  experimental  data,  however,  in  this 
case  Vs(t)  can  only  be  calculated  by  dividing  the  distance  between  two  adjacent  transducers  by 
the  time  it  takes  for  the  shock  to  cross  both  transducers.  Note  that  this  is  equivalent  to  differ¬ 


entiating  Xs(t).  The  corresponding  velocity  history  is  shown  in  Fig.  2i. 

In  the  authors'  opinion,  Xs(t)  and  Vs(t)  are  the  most  comprehensive  and  accurate  char¬ 
acterization  of  the  unsteadiness  in  these  interactions.  Once  they  are  determined,  Xs(t)  and  Vs(t) 
are  analyzed  using  standard  time-series  analysis  techniques  and  are  correlated  with  wall  pres¬ 
sure  measurements  made  upstream  and  downstream  of  the  intermittent  region  or  with  other 
conditionally  sampled  signals  obtained  from  wall  pressure  measurements  in  a  given  interaction. 
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3.3.  Derivation  of  instantaneous  pressure  ratio  history 

In  addition  to  Xs(t)  and  Vs(t),  it  is  also  desirable  to  know  the  histories  of  the  ratio  of 
static  pressures  across  the  separation  shock.  To  obtain  this,  one  would  ideally  need  an  instru¬ 
ment  which  traveled  with  the  separation  shock.  Clearly,  it  is  not  possible  to  make  these  inea- 
surements.  However,  an  approximation  of  the  static  pressure  ratio  history  can  be  obtained 
from  simultaneous  measurements  of  wall  pressure  fluctuations  in  the  intermittent  region. 

To  calculate  the  instantaneous  pressure  ratio  across  the  separation  shock,  Rs(t),  one 
must  first  determine  its  position  at  that  instant  in  time.  Note  that  the  separation  shock  position 
history,  Xs(t),  has  already  been  determined  and  is  shown  in  Fig.  2h.  At  any  time,  Rs(t)  is 
given  by  P2(t)/Pi(t),  where  P2(t)  are  the  pressure  fluctuations  measured  by  the  transducer  im¬ 
mediately  downstream  of  the  separation  shock  and  Pi(t)  are  those  fluctuations  measured  by  the 
upstream  transducer.  As  can  be  seen  in  Fig.  2h,  the  separation  shock  is  in  bin  5  from  time 
t  =  rs  to  t  =  r4  .  During  these  times,  Pi(t)  and  P2(t)  fluctuations  are  measured  by  channels 
4  and  5,  respectively,  and  Rs(t)  is  calculated  from  these  pressure  fluctuations.  Similarly,  for 
r3  <  t  <  fs  the  separation  shock  is  in  bin  3,  and  Rs(t)  is  calculated  from  pressure  fluctuations 
measured  by  channels  2  and  3.  As  a  result  of  zero-drift  problems,  boundary  layer  pressure 
level,  Pwo,  will  differ  from  channel  to  channel.  This  can  lead  to  inaccurate  calculation  of  Rs(t). 
Therkore,  to  avoid  any  bias  in  Rs(t),  measured  Pwo  are  first  subtracted  from  upstream,  Pi(t) 
and  downstream,  P2(t)  pressure  fluctuations  and  freestream  static  pressure,  Poo,  is  added  to 
both.  Corresponding  pressure  ratio  fluctuations,  Rs(t),  calculated  from  data  shown  in  Figs.  2a 

through  2f,  are  plotted  in  Fig.  2j.  r,  /  ^  • 

It  should  be  noted  that,  when  the  separation  shock  is  in  bin  1,  Pi(t)  fluctuations  are  not 
known,  and  Rs(t)  cannot  be  calculated.  Similarly,  when  the  shock  is  in  bin  9,  P2(t)  fluctua¬ 
tions  are  not  known  (see  the  schematic).  Those  times  during  which  Rs(t)  cannot  be  calculated 
are  ignored  in  the  analysis. 


3.4.  Variable  interval  cross-correlation 

Cross-correlation  of  two  simultaneously  sampled  signals  is  generally  calculated  using 
standard  time-series  analysis  techniques.  The  data  are  first  divided  into  equal  interval  records, 
which  are  Fourier  transformed  independently.  The  transformed  results  are  then  ensemble- 
averaged  to  give  cross-spectral  density  function.  The  cross-correlation  result  is  obtained  by 
taking  the  inverse  Fourier  transform  of  the  cross-spectral  density  function. 

In  some  cases,  particularly  when  conditional  sampling  analysis  is  used,  it  is  necessary 
to  work  with  variable  interval  data  blocks.  The  Variable  Interval  Cross-Correlation  algorithm 
has  been  developed  for  exactly  this  purpose.  This  algorithm  performs  an  algebraic  calculation 
of  cross-correlation  between  two  simultaneously  sampled  channels  of  data.  In  this  rriethod,  the 
data  from  the  first  channel  are  divided  into  a  certain  number  of  blocks,  each  of  which  can  be 
different  in  size.  Each  data  block  is  then  correlated  with  the  corresponding  data  block  extracted 
from  the  second  channel.  This  is  repeated  for  all  data  blocks  and  the  final  cross-correlation 
result  is  obtained  by  ensemble-averaging  the  individual  cross-correlation  results  for  each  data 
block.  To  improve  accuracy,  the  data  block  from  the  second  channel  is  padded  on  both  sides. 

The  results  presented  in  this  paper  have  been  calculated  using  both  the  standard  and  the 
Variable  Interval  Cross-Correlation  methods. 


3.5.  Digital  filter 

A  Finite  Impulse  Response,  FIR,  digital  filter  has  been  used  to  examine  band-limited 
correlations  between  separation  shock  dynamics  (i.e.,  Xs(t)  and  Vs(t))  and  other  measure¬ 
ments.  To  do  this,  the  signals  were  digitally  filtered  in  either  the  low-pass  or  band-pass  modes 
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prior  to  any  analysis.  Correlations  of  Xs(t)  and  Vs(t)  with  fluctuating  pressure  measurements 
were  examined  in  frequency  bands  of  0-50  kHz,  0-3  kHz  and  3-50  kHz.  It  should  be  noted 
that,  Vs(t)  was  itself  never  filtered,  instead  the  so-called  filtered  velocity  fluctuations,  Vs(t) 
were  merely  the  derivative  of  the  filtered  Xs(t)  signal.  A  piecewise  smooth  signal,  such  as  , 
can  be  represented  by  the  sum  of  an  infinite  series  of  sinusoids.  When  the  signal  is  low-pass 
filtered,  the  result  corresponds  essentially  to  a  truncated  series  of  sinusoids.  Consequently,  the 
filtered  signal  displays  the  well  known  Gibbs  phenomenon,  however,  the  effect  of  these 
“ripples”  is  negligible  for  the  types  of  analysis  performed  in  this  study. 

3.6.  CFD  model 

In  addition  to  the  analysis  of  the  experimental  data,  some  of  the  physical  phenomena 
have  been  examined  using  an  inviscid,  quasi  one-dimensional  computational  model.  In  this 
model,  time-accurate  one-dimensional  Euler  equations  are  used  to  solve  supersonic/subsonic 
inviscid  flow  in  a  diverging  channel.  The  inlet  and  outlet  boundary  condition  were  so  pre¬ 
scribed  to  force  a  normal  shock  to  occur  in  the  mid-section  of  the  channel.  Once  the  steady- 
state  solution  was  reached,  the  inlet  and/or  outlet  boundary  conditions  were  varied  using  sinu¬ 
soidal  perturbations  whose  amplitude  and  frequency  could  be  controlled  independently.  The 
response  of  the  normal  shock  to  these  perturbations  was  monitored  by  recording  its  position  at 
each  time  step.  Although  this  numerical  model  does  not  encompass  all  physical  phenomena 
observed  in  experiments,  it  does  nevertheless  model  the  unsteadiness  of  a  shock  and  convec¬ 
tion  of  disturbances  in  an  inviscid  flow.  It  can,  therefore,  be  used  to  show  qualitatively  how  a 
shock  will  respond  to  perturbations  in  the  boundary  conditions. 

It  should  be  emphasized  that,  the  authors  do  not  claim  to  have  modeled  the  physical 
phenomena  observed  in  experiments.  Instead,  this  CFD  model  has  been  introduced  as  a  tool  to 
examine  some  aspects  of  the  unsteadiness  observed  in  shock  wave  /  turbulent  boundary  layer 
interactions.  It  will  be  shown  later  that  the  results  obtained  from  this  numerical  analysis  sup¬ 
port  one  part  of  the  mechanism  of  unsteadiness  deduced  from  the  analysis  of  the  experimental 
data. 

4.  Discussion  of  Results 

To  recap,  the  objective  of  the  current  work  is  to  examine  the  separation  shock  dynamics 
in  several  interactions  and  determine  whether  a  universal  model  describing  the  unsteadiness  can 
be  developed.  To  do  this,  the  unsteadiness  is  characterized  in  terms  of  separation  shock  posi¬ 
tion,  Xs(t)  and  velocity,  Vs(t)  histories.  These  quantities  have  been  correlated  with  condition¬ 
ally  sampled  pressure-ratio  fluctuations  and  wall  pressure  fluctuations  measured  upstream 
and/or  downstream  of  the  intermittent  region.  Before  discussing  the  results  of  these  correla¬ 
tions,  some  basic  measurements  describing  these  interactions  will  be  presented  first. 
Following  the  discussion  of  cross-correlation  results,  a  model  describing  the  physical  causes  of 
separation  shock  unsteadiness  will  be  developed. 

4.1.  Mean,  standard  deviation,  and  intemiittency  distributions 

Distributions  of  normalized  mean  wall  pressure,  Pw  /  Poc  are  plotted  versus  the 
normalized  distance,  x/5o,  in  Fig.  3a.  It  should  be  noted  that  6o  is  not  the  appropriate  normal- 
izer  in  all  cases.  In  fact,  the  data  from  the  25  deg  swept  compression  ramp  interaction  should 
be  plotted  in  the  conical  coordinate  system  (Erengil  &  Dolling  1993),  and  those  from  the  blunt- 
fin  interactions  normalized  by  leading  edge  diameter.  However,  the  intent  here  is  to  demon¬ 
strate  the  similarities  and  differences  of  the  distributions  from  the  unswept  compression  ramp 
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(UCR),  the  25  deg  swept  compression  ramp  (SCR-25),  the  8  deg  swept  blunt-fin  (SBF-8), 
and  the  30  deg  swept  blunt-fin  (SBF-30)  interactions.  The  separation  lines,  S,  and  upstream 
influence  lines,  UI,  for  the  UCR  and  SBF-8  interactions  are  shown  in  Fig.  3a,  and  those  for 
SCR-25  and  SBF-30  are  shown  in  Fig.  3b.  As  expected,  Pw  /  distributions  are  similar  in 
that  they  all  increase  gradually  through  the  intermittent  regions  upstream  of  the  compression 
comer  or  the  blunt-fin  root.  Note  that  the  spatial  extent  of  the  SBF-30  interaction  is  about  one- 
third  of  the  SBF-8  interaction.  Similarly,  the  upstream  influence  of  the  SCR-25  interaction  is 
smaller  than  that  of  the  UCR  interaction. 

The  normalized  standard  deviation,  Opw/Poo,  distributions  are  plotted  in  Fig.  3b.  Since 
the  local  maxima  in  these  distributions  are  generated  by  the  motion  of  the  separation  shock  in 
the  intermittent  region,  their  magnitudes  provide  a  measure  of  the  strength  of  the  separation 
shock.  Based  on  the  relative  magnitudes  of  local  maxima  in  apw/P  oo  shown  in  Fig.  3b,  the 
separation  shock  has  almost  the  same  strength  in  the  UCR  and  SBF-8  interactions  whereas  it  is 
slightly  weaker  in  the  SBF-30  and  SCR-25  interactions. 

Finally,  the  intermittency  distributions  for  these  interactions  are  plotted  in  the  same 
coordinates  in  Fig.  3c.  Since  the  intermittency,  y,  has  a  value  between  0  and  1,  the  slope  of 
these  curves  provides  a  simple  measure  of  the  length  of  the  intermittent  region.  By  visual 
inspection  of  the  results  in  Fig.  3c,  it  can  easily  be  concluded  that  the  intermittent  regions  in  the 
UCR  and  SBF-8  interactions  are  about  the  same  length.  The  SBF-30  interaction  has  a  smaller 
intermittent  region  and  that  in  the  SCR-25  interaction  is  the  smallest. 

4.2.  Power  spectral  densities  ofXs(t)  and  Vs(t) 

Power  spectral  density  coefficients  of  Xs(t)  and  Vs(t)  are  plotted  in  normalized  form  in 
Figs.  4a  and  4b,  respectively.  These  specific  cases  are  identified  in  the  legend  in  Fig.  4a,  and 
for  clarity,  those  results  in  Fig.  4b  have  been  offset  by  the  amounts  indicated  in  its  legend. 
When  plotted  in  this  form,  the  total  area  under  the  normalized  distributions  is  unity,  and  the 
area  within  a  frequency  band  represents  the  contribution  to  the  variance  from  that  frequency 
band. 

As  can  be  seen  in  Fig.  4a,  there  is  practically  no  contribution  to  the  variance  of  Xs(t) 
from  frequencies  in  excess  of  about  10  kHz  in  any  of  these  interactions.  To  verify  that  the 
transducer  spacing  does  not  impose  a  spatial  filter  on  the  Xs(t)  signal,  the  authors  used  half  as 
many  transducers  at  double  spacing  to  characterize  the  separation  shock  position  history. 
There  were  few  noticeable  differences  in  the  power  spectral  density  distributions. 
Consequently,  changes  in  the  separation  shock  position  occur  at  frequencies  below  10  kHz  in 
all  these  interactions.  Clearly,  dominant  frequencies  in  Xs(t)  increase  progressively  from  about 
0.2  kHz  in  the  UCR  interaction  to  about  3  kHz  in  the  SCR-25  interaction.  In  the  blunt-fin 
interactions,  dominant  frequencies  in  Xs(t)  are  between  these  two  extremes  and  are  higher  in 
the  SBF-30  interaction  than  they  are  in  the  SBF-8  interaction.  Increasing  dominant  frequencies 
are  consistent  with  decreasing  intermittent  region  lengths  (Gonsalez  &  Dolling  1993). 

Several  observations  can  be  made  from  normalized  power  spectral  density  distributions 
shown  in  Fig.  4b.  First,  the  sharp  drop  at  50  kHz  is  due  to  digital  filtering  of  the  separation 
shock  position  history,  Xs(t).  As  was  mentioned  earlier,  Vs(t)  is  never  filtered,  and  the  so- 
called  filtered  velocity  fluctuations  are  simply  the  derivative  of  the  filtered  Xs(t)  signal.  The 
effect  of  the  digital  filter  could  not  be  seen  in  Fig.  4a,  because  there  is  practically  no  energy  in 
Xs(t)  at  frequencies  above  50  kHz.  Second,  the  frequency  content  of  the  velocity  fluctuations 
is  essentially  broadband  particularly  in  the  UCR  and  SBF-8  interactions.  In  the  SBF-30  inter¬ 
action,  the  frequency  distribution  is  bimodal  with  a  dominant  peak  in  the  3  to  4  kHz  range.  In 
the  SCR-25  interaction,  dominant  frequencies  are  within  the  range  of  4  to  5  kHz  and  very  well 
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pronounced.  Third,  there  is  a  gradual  and  a  subtle  trend  towards  higher  frequencies  in  Vs(t) 
with  decreasing  intermittent  region  length. 

To  summarize,  the  results  presented  so  far  show  the  similarities  and  differences  in  all 
the  interactions  studied.  As  indicated  by  the  mean  results  presented  in  Fig.  3,  although  their 
scales  differ,  these  interactions  all  have  a  similar  separation  process.  However,  dominant  fre¬ 
quencies  characterizing  the  separation  shock  unsteadiness  are  quite  different  in  each  of  these 
interactions. 


4.3.  Correlation  ofVs(t)  with  pressure  ratio  fluctuations 

For  an  unsteady  shock  wave,  the  governing  equations  can  be  reduced  to  an  expression 
in  which  the  shock  speed  is  given  by  the  ratio  of  static  quantities  across  the  shock,  primarily 
the  static  pressure  ratio.  Hence,  it  is  useful  to  examine  whether  there  is  a  correlation  between 
separation  shock  velocity,  Vs(t)  and  conditionally  sampled  wall  pressure  ratio,  Rs(t)  histories. 
As  was  discussed  earlier,  Rs(t)  can  be  derived  from  simultaneous  measurements  of  multi¬ 
channel  wall  pressure  fluctuations  made  in  the  intermittent  region. 

Normalized  results  of  the  correlation  of  Vs(t)  and  Rs(t)  are  shown  in  Fig.  5.  For  clar¬ 
ity,  results  from  the  SBF-8,  SBF-30,  and  SCR-25  interactions  have  been  offset  as  indicated  in 
the  figure.  The  negative  correlation  peak  implies  that  negative  fluctuations  in  separation  shock 
velocity,  Vs(t)  correlate  with  positive  fluctuations  in  wall  pressure  ratio,  Rs(t)  and  vice  versa. 
Furthermore,  a  negative  time  delay  means  that  fluctuations  in  Vs(t)  are  preceded  by  those  in 
Rs(t).  Thus  an  increase  in  the  static  pressure  ratio  across  the  separation  shock  causes  a 
decrease  in  its  velocity  and  vice  versa.  It  should  be  noted  that  a  change  in  separation  shock 
velocity  does  not  necessarily  imply  a  change  in  its  direction  of  motion  —  it  simply  means  that 
the  separation  shock  is  either  speeding  up  or  slowing  down.  The  direction  of  motion  of  the 
separation  shock  is  given  by  the  sign  of  its  velocity,  and  in  this  case  a  negative  velocity  means 
that  the  separation  shock  is  moving  in  the  upstream  direction  and  vice  versa.  The  mechanism 
of  unsteadiness,  therefore,  is  an  integral  result  of  fluctuations  in  separation  shock  velocity, 
where  positive  shock  velocities  integrate  to  downstream  excursions  and  vice  versa.  Finally, 
these  results  show  that  this  mechanism  is  common  to  and  identical  in  all  interactions. 

To  the  best  of  the  authors'  knowledge,  these  results  provide  the  first  experimental  evi¬ 
dence  of  a  cause-and-effect  relationship  between  wall  pressure  fluctuations  and  the  separation 
shock  wave  unsteadiness  in  these  interactions.  To  reiterate,  fluctuations  in  the  static  pressure 
ratio  across  the  separation  shock  induce  fluctuations  in  its  velocity  which  in  turn  cause  a  change 
in  separation  shock  position.  It  should  be  noted,  however,  that  the  separation  shock  unsteadi¬ 
ness  described  by  this  mechanism  is  the  small-scale,  jittery  motion,  since  Rs(t)  is  obtained  by 
“following”  the  separation  shock.  In  other  words,  since  this  analysis  focuses  on  those  fluctua¬ 
tions  in  the  vicinity  of  the  separation  shock  foot,  it  cannot  account  for  the  large-scale,  global 
motion  of  the  separation  shock.  To  determine  the  cause  of  the  large-scale  unsteadiness,  sepa¬ 
ration  shock  dynamics  must  be  correlated  with  upstream  and  downstream  pressure  fluctuations 
separately.  These  results  and  others  are  presented  in  the  following  sections. 

4.4.  Correlation  ofVs(t)  with  upstream  pressure  fluctuations 

Normalized  results  of  the  correlation  of  Vs(t)  with  upstream  wall  pressure  fluctuations 
are  shown  in  Fig.  6.  For  clarity,  results  from  the  SBF-8,  SBF-30  and  SCR-25  interactions 
have  been  offset  as  indicated  in  the  figure.  Note  that  the  normalized  coefficients  for  these  cor¬ 
relations  are  much  lower  than  those  shown  in  Fig.  5.  This  is  because  the  separation  shock 
velocity  fluctuations  actually  correlate  with  fluctuations  in  Rs(t),  and  these  may  or  may  not  be 
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affected  by  the  upstream  pressure  fluctuations.  Obviously,  the  static  pressure  ratio  is  also  a 
function  of  the  instantaneous  value  of  the  downstream  pressure  fluctuations. 

Regardless,  there  is  a  correlation  and  the  results  show  a  positive  peak  at  negative  time 
delay.  This  suggests  that  a  positive  (negative)  fluctuation  in  the  upstream  pressure  can  gener¬ 
ate  a  negative  (positive)  fluctuation  in  static  pressure  ratio,  Rs(t)  which  causes  a  positive 
(negative)  fluctuation  in  Vs(t).  The  peak  in  the  SCR-25  interaction  is  the  smallest,  but  this 
should  be  expected  since  this  interaction  is  highly  swept  and  some  turbulent  structures  which 
cross  over  the  upstream  transducer  will  not  reach  the  separation  shock  foot  along  the  line  of 
measurement  which  is  perpendicular  to  the  loeal  separation  line.  These  results  are  consistent 
with  and  support  those  presented  in  Fig.  5.  In  each  of  the  blunt-fin  interactions,  there  is  also  a 
secondary  positive  peak  at  an  earlier  time-delay.  This  is  due  to  the  strong  recirculation  present 
in  blunt-fin  interactions  (Kleifges  &  Dolling  1993).  Some  turbulent  structures  which  have 
crossed  over  the  upstream  transducer  travel  to  the  fin  root  and  back,  and  generate  a  fluctuation 
in  Rs(t)  which  in  turn  causes  a  fluctuation  in  Vs(t). 

To  demonstrate  this  physical  phenomenon  further,  wall  pressure  measurements  made 
upstream  and  downstream  of  the  intermittent  regions  of  these  interactions  have  been  correlated 
with  each  other.  Normalized  results  of  these  correlations  are  shown  in  Fig.  7.  For  clarity,  the 
results  from  the  SBF-8,  SBF-30  and  SCR-25  interactions  have  been  offset  as  indicated  in  the 
figure.  It  should  be  noted  that  these  results  are  obtained  from  actual  wall  pressure  measure¬ 
ments  and  not  from  conditionally  generated  signals.  In  the  UCR  interaction  there  is  only  a 
single  peak  at  positive  time  delay.  This  is  because  the  turbulent  stmctures  which  cross  over  the 
upstream  transducer  are  convected  through  the  intermittent  region  and  also  cross  over  the 
downstream  transducer,  but  are  not  recirculated.  A  similar  convective  process  is  also  evident 
in  the  SCR-25  interaction.  Furthermore,  the  lower  correlation  magnitude  suggests  that  some  of 
the  turbulent  structures  which  have  crossed  over  the  upstream  transducer  are  swept  away 
before  they  reach  the  downstream  transducer.  It  should  be  noted  that  differences  in  the  timing 
of  the  primary  peaks  are  a  direct  result  of  different  spacings  between  the  two  transducers.  In 
the  SBF-8  and  SBF-30  interactions  there  is  a  primary  and  a  secondary  peak  at  positive  time 
delays.  The  primary  peak  is  due  to  the  turbulent  structures  which  have  crossed  over  the 
upstream  transducer,  convected  through  the  intermittent  region  and  also  crossed  over  the 
downstream  transducer.  The  secondary  peak,  on  the  other  hand,  is  caused  by  those  structures 
which  are  recirculated  within  the  separation  zone,  because  the  time  delay  for  each  of  these  sec¬ 
ondary  peaks  corresponds  approximately  to  the  time  required  for  turbulent  structures  to  travel 
to  the  root  of  the  blunt-fin  and  back  over  the  downstream  transducer.  In  addition  to  spatial 
decay  of  turbulent  structures,  the  lower  correlation  magnitude  for  the  secondary  peaks  also 
suggests  that  only  some  of  the  structures  are  recirculated.  Furthermore,  the  time  difference 
between  the  two  peaks  indicates  that  the  recirculating  structures  require  longer  time  in  the  SBF- 
8  interaction  than  they  do  in  the  SBF-30  interaction.  This  is  because  the  latter  is  a  much 
smaller  interaction  than  the  former,  as  can  be  ascertained  from  Fig.  3. 

To  recap,  the  results  presented  in  Fig.  6  indicate  that  correlation  of  velocity  fluctuations 
with  upstream  pressure  fluctuations  are  essentially  the  same  for  all  the  interactions  studied. 
The  secondary  peaks  in  blunt-fin  interactions  reflect  some  characteristics  unique  to  these  inter¬ 
actions  as  explained  above.  Furthermore,  although  these  results  provide  additional  support  for 
the  mechanism  developed  from  correlations  presented  in  Fig.  5,  they  clearly  do  not  correlate 
with  the  large-scale,  global  motion  of  the  separation  shock.  In  other  words,  turbulent  fluctua¬ 
tions  in  the  incoming  boundary  layer  can  only  induce  small-scale  unsteadiness  by  generating 
fluctuations  in  the  ratio  of  static  quantities  across  the  separation  shock.  Consistent  with  the 
approach  outlined  at  the  end  of  the  previous  section,  results  from  correlation  of  separation 
shock  dynamics  and  downstream  pressure  fluctuations,  as  well  as  others,  are  presented  in  the 
following  section. 
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4.5.  Correlation  ofVs(t)  with  downstream  pressure  fluctuations 

The  interpretation  of  the  correlation  of  separation  shock  velocity,  Vs(t)  and  downstream 
pressure  fluctuations  is  somewhat  complicated,  since  these  fluctuations  are  coupled. 
Normalized  results  of  the  correlation  of  Vs(t)  and  wall  pressure  fluctuations  are  plotted  in  Fig. 
8  for  all  interactions  studied.  For  clarity,  the  results  from  the  SBF-8,  SBF-SO  and  SCR-25 
interactions  have  been  offset  as  indicated  in  the  figure.  In  all  four  cases,  the  downstream 
transducer  is  located  at  or  just  downstream  of  the  separation  line.  Clearly,  the  results  are  dif¬ 
ferent  for  different  interactions.  For  the  SCR-25  interaction,  for  example,  there  is  a  well 
defined  positive  peak  at  negative  time-delay  and  a  negative  peak  at  positive  time-delay.  This 
means  that,  positive  (negative)  fluctuations  in  pressure  precede  positive  (negative)  fluctuations 
in  velocity,  and  that  negative  (positive)  fluctuations  in  velocity  are  followed  by  positive 
(negative)  fluctuations  in  pressure.  This  implies  that  there  is  a  cyclic  relationship  between  sep¬ 
aration  shock  velocity  fluctuations  and  wall  pressure  measurements  near  the  separation  line. 

Similar  observations  are  also  valid  for  the  SBF-30,  SBF-8,  and  UCR  interactions, 
although  the  correlation  peaks  broaden  and  shift  to  consistently  greater  time  delays.  This  im¬ 
plies  that  characteristic  (or  dominant)  frequencies  of  this  cyclic  relationship  between  separation 
shock  velocity  and  downstream  pressure  fluctuations  shift  progressively  to  lower  frequency 
bands  for  the  SBF-30,  SBF-8,  and  UCR  interactions,  respectively.  It  should  be  noted  that 
these  trends  are  consistent  with  those  observed  in  power  spectral  density  coefficients  of  sepa¬ 
ration  shock  position,  Xs(t)  and  velocity,  Vs(t)  histories.  Furthermore,  as  the  dominant  fre¬ 
quencies  of  this  cyclic  relationship  decrease,  other  correlation  peaks  around  'T=0  become 
apparent  in  at  least  the  UCR  and  SBF-8  interactions.  These  observations  suggest  that  there  are 
two  different  physical  phenomena  by  which  velocity  fluctuations  correlate  with  downstream 
wall  pressure  fluctuations,  and  that  these  phenomena  have  different  characteristic  frequencies. 
Therefore,  to  better  understand  the  results  of  correlation  between  these  fluctuations,  the  signals 
have  been  digitally  filtered  in  low-pass  and  high-pass  modes  to  examine  contributions  from 
different  frequency  bands. 

Normalized  correlation  of  band-pass  filtered  fluctuations  are  shown  in  Fig.  9  for  these 
interactions.  The  cut-off  frequencies  were  set  at  3  and  50  kHz.  Again,  results  from  the  SBF- 
8,  SBF-30  and  SCR-25  interactions  have  been  offset  for  clarity.  In  this  frequency  band,  the 
correlation  results  are  qualitatively  the  same  for  all  the  interactions  studied.  The  dominant  peak 
in  these  correlations  has  a  negative  magnitude  and  occurs  at  negative  time  delay.  This  suggests 
that  negative  (positive)  fluctuations  in  downstream  pressure  precede  positive  (negative)  fluc¬ 
tuations  in  separation  shock  velocity.  These  results  are  entirely  consistent  with  and  provide 
additional  support  for  the  mechanism  described  in  the  discussion  of  correlations  presented  in 
Fig.  5.  Note  that  a  negative  fluctuation  in  downstream  pressure  may  generate  a  negative  fluc¬ 
tuation  in  Rs(t)  which  causes  a  positive  fluctuation  in  separation  shock  velocity.  Similarly,  a 
positive  fluctuation  in  downstream  pressure  may  increase  the  pressure  ratio  across  the  separa¬ 
tion  shock  and  thereby  induce  a  negative  velocity.  The  magnitude  of  these  correlatiori  peaks  is 
smaller  than  those  shown  in  Fig.  5,  because  not  all  positive  (negative)  fluctuations  in  down¬ 
stream  pressure  will  generate  a  positive  (negative)  fluctuation  in  Rs(t).  The  other  peaks  in  the 
correlation  results  are  due  to  the  strong  coupling  between  fluctuations  in  downstream  pressure 
and  those  in  separation  shock  velocity  as  well  as  in  position  histories.  It  should  be  noted  that 
separation  shock  strength  increases  as  the  shock  moves  downstream,  and  vice  versa  (Erengil  & 
Dolling  1991a).  Therefore,  a  positive  (negative)  fluctuation  in  separation  shock  velocity, 
which  causes  a  downstream  (upstream)  motion,  results  in  a  stronger  (weaker)  shock,  which  in 
turn  generates  a  positive  (negative)  fluctuation  in  downstream  pressure.  Clearly,  the  strong 
coupling  between  separation  shock  dynamics  and  downstream  fluctuations  make  the  interpreta¬ 
tion  of  these  results  quite  difficult. 
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To  further  examine  flow  dynamics  in  the  downstream  regions  of  these  interactions, 
wall  pressure  measurements  under  the  separated  bubble  have  been  correlated  with  each  other. 
Normalized  results  of  these  correlations  are  shown  in  Fig.  10.  Again,  these  results  have  been 
offset  for  clarity  as  indicated  in  the  figure.  To  minimize  the  effect  of  shock  unsteadiness  on 
these  results,  care  was  taken  in  the  analysis  to  ensure  that  no  data  were  included  in  which  the 
separation  shock  crossed  over  any  of  these  transducers.  Furthermore,  the  pressure  signals 
were  digitally  filtered  with  the  cut-off  frequencies  set  at  3  and  50  kHz.  This  analysis,  there¬ 
fore,  focuses  on  the  turbulent  fluctuations  and  the  results  represent  the  dynamics  of  turbulent 
stmctures.  In  the  UCR  interaction  there  is  a  positive  peak  occurring  at  positive  time  delay. 
This  suggests  that  the  turbulent  structures  are  convected  in  the  downstream  direction,  and  if 
there  is  backflow,  it  is  not  strong  enough  to  be  detected  by  wall  pressure  measurements.  In  the 
SBF-8,  SBF-30  and  SCR-25  interactions,  there  are  peaks  at  positive  and  negative  time  delays, 
indicating  that  turbulent  structures  are  convected  in  both  upstream  and  downstream  directions 
in  the  separated  region.  These  observations  further  support  the  correlation  results  presented  in 
Fig.  9,  and  show  clearly  that  turbulent  structures  can  be  convected  upstream  to  the  vicinity  of 
the  separation  shock,  change  the  ratio  of  the  static  quantities  and  thereby  induce  a  fluctuation  in 
the  shock  velocity,  Vs(t). 

Normalized  correlations  of  low-pass  filtered  fluctuations  are  shown  in  Fig.  1 1  for  the 
UCR,  SBF-8,  SBF-30  and  SCR-25  interactions.  These  results  correspond  to  those  presented 
in  Fig.  8  except  that  they  have  been  low-pass  filtered  at  3  kHz.  These  results  have  also  been 
offset  as  indicated  in  the  figure.  In  this  frequency  band,  the  results  are  qualitatively  the  same 
for  all  the  interactions.  They  all  have  a  positive  peak  at  a  negative  time  delay  and  a  negative 
peak  at  a  positive  time  delay.  The  difference  is  in  the  broadening  of  the  peaks  which  suggests  a 
progressive  shift  in  dominant  frequencies  towards  the  lower  end  of  the  0-3  kHz  range  with  the 
lowest  frequencies  present  in  the  UCR  interaction.  The  peaks  at  negative  time  delay  suggest 
that  positive  (negative)  fluctuations  in  downstream  pressure  precede  positive  (negative)  fluc¬ 
tuations  in  separation  shock  velocity,  whereas  those  at  positive  time  delay  suggest  that  positive 
(negative)  fluctuations  in  pressure negative  (positive)  fluctuations  in  velocity.  Before 
describing  this  correlation,  it  will  be  helpful  to  note  that  the  wall  pressure  levels  at  these  down¬ 
stream  stations  are  higher  when  the  separation  shock  is  at  the  upstream  end  of  the  interniittent 
region  and  lower  when  it  is  at  the  downstream  end  (see  also  Fig.  12).  Thus  a  shock  motion  in 
the  upstream  direction  (i.e.,  a  negative  velocity)  will  be  followed  by  an  increase  in  pressure  at 
the  downstream  station  and  vice  versa.  This  is  entirely  consistent  with  negative  peaks  at  posi¬ 
tive  time  delays  seen  in  Fig.  11.  That  there  exist  positive  peaks  at  negative  time  delays  suggest 
that  this  correlation  is  a  cyclic  process.  In  other  words,  a  high  downstream  pressure  level  is 
followed  by  a  shock  motion  in  the  downstream  direction,  which  leads  to  a  decrease  in  the 
downstream  pressure  level,  which  in  turn  is  followed  by  a  shock  motion  in  the  upstream  direc¬ 
tion,  leading  to  an  increase  in  the  downstream  pressure  level.  Clearly,  this  is  a  cyclic  process 
which  makes  it  difficult  to  establish  a  cause-and-effect  relationship.  However,  it  does  suggest 
that  there  is  a  correlation  between  large-scale  separation  shock  motion  and  expansion  and  con¬ 
traction  of  the  separation  bubble,  since  an  upstream  motion  of  the  shock  corresponds  to 
expansion  of  the  bubble  whereas  a  downstream  motion  corresponds  to  contraction  of  the 
bubble. 


4.6.  Correlation  ofXs(t)  with  downstream  pressure  fluctuations 

The  mechanism  suggested  by  the  correlation  between  separation  shock  motion  and 
expansion  and  contraction  of  the  separation  bubble  is  quite  different  from  that  dictated  by  the 
ratio  of  static  quantities  across  the  separation  shock  foot.  In  the  latter,  the  shock  motion  is  an 
integral  result  of  velocity  fluctuations  caused  by  perturbations  in  the  ratio  of  static  quantities 
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across  the  shock.  In  the  former,  the  separation  shock,  which  forms  inunediately  upstream  of 
the  separation  bubble,  is  simply  displaced  by  such  physical  processes  as  the  expansion  and 
contraction  of  the  bubble.  Thus  in  the  former  case,  the  relevant  parameter  characterizing  the 
separation  shock  unsteadiness  is  its  position  and  not  velocity  history.  Therefore,  it  is  Xs(t) 
which  should  be  correlated  with  measurements  quantifying  the  dynamics  of  expansion  and 
contraction  of  the  separation  bubble.  Following  is  a  short  discussion  explaining  how  the 
expansion  and  contraction  of  the  bubble  can  best  be  quantified  by  wall  pressure  measurements. 

Figure  12  (Erengil  &  Dolling  1991a)  shows  ensemble-averaged  pressure  distributions 
for  different  separation  shock  positions  in  the  intermittent  region  of  the  UCR  interaction.  The 
solid  line  is  the  overall  time-averaged  pressure  distribution,  and  the  dashed  lines  correspond  to 
different  shock  positions  as  shown  in  the  figure.  The  most  upstream  dashed  line  distribution, 
indicated  by  open  square  markers,  corresponds  to  a  large  separation  bubble,  whereas  that, 
indicated  by  filled  diamond  markers,  corresponds  to  a  small  separation  bubble.  These  distri¬ 
butions  show  that  expansion  and  contraction  of  the  separation  bubble  can  be  quantified  by  wall 
pressure  measurements  made  in  the  vicinity  of  the  separation  line.  It  should  he  noted  that  this 
is  the  optimum  location  for  monitoring  expansion  and  contraction  of  the  separation  bubble, 
since  measurements  made  under  the  separation  bubble,  say  at  x/5o=-0.5  (see  Fig.  12),  are 
essentially  independent  of  the  size  of  the  bubble,  and  those,  say  at  40%  intermittency,  cannot 
quantify  the  dynamics  of  the  separation  for  60%  of  the  time  because  the  transducer  will  be 

under  the  undisturbed  boundary  layer  for  these  times. 

Figure  13  shows  normalized  results  of  correlation  between  separation  shock  position, 
Xs(t)  and  wall  pressure  measurements  near  the  separation  line  in  the  UCR,  SBF-8,  SBF-30 
and  SCR-25  interactions.  For  clarity,  results  from  the  SBF-8,  SBF-30  and  SCR-25  interac¬ 
tions  have  been  offset  as  indicated  in  the  figure.  Note  that  these  results  represent  correlations 
in  the  entire  frequency  range,  that  is  0  to  50  kHz.  First,  the  results  are  qualitatively  the  same  in 
all  the  interactions.  The  difference  is  in  the  broadening  of  the  peaks  which  is  also  accompanied 
by  a  decrease  in  their  magnitudes.  As  was  discussed  earlier,  the  broadening  of  the  peaks  indi¬ 
cates  a  decrease  in  the  characteristic  frequencies  of  the  unsteadiness  which  is  also  ch^acterized 
by  power  spectral  density  distributions  of  separation  shock  position  histories  shown  in  Fig.  4a. 
Most  importantly,  these  correlations  have  a  strong  negative  peak  at  negative  time  delay.  This 
means  that  negative  (positive)  fluctuations  in  pressure  precede  positive  (negative)  fluctuations 
in  separation  shock  position.  It  should  be  noted  that  higher  levels  of  pressure  near  the  separa¬ 
tion  line  correspond  to  a  larger  separation  bubble  and  vice  versa  (also  see  Fig.  12).  Given  this, 
the  results  presented  in  Fig.  13  show  that  the  expansion  of  the  separation  bubble  causes  a 
shock  motion  in  the  upstream  direction,  whereas  the  contraction  of  the  bubble  results  in  a 
motion  in  the  downstream  direction. 


4. 7.  Physical  model  of  the  unsteadiness 

Based  on  the  experimental  results  presented  in  previous  sections  a  physical  model 
describing  the  separation  shock  unsteadiness  in  all  four  shock  wave  /  turbulent  boundary  layer 
interactions  can  be  constructed.  In  this  model,  separation  shock  motion,  can  be  caused  by  two 
different  physical  phenomena  which  comprise  the  two  parts  of  the  mechanism.  In  the  first 
part,  the  unsteadiness  is  caused  by  perturbations  in  the  ratio  of  static  quantities  across  the  sepa¬ 
ration  shock.  The  ratio  of  static  pressure,  P2(t)/Pi(t),  for  example,  can  be  changed  by 
fluctuations  in  P2(t),  Pi(t)  or  both.  This,  according  to  governing  equations,  will  cause  either 
an  increase  or  a  decrease  in  separation  shock  velocity,  which  in  turn  integrates  to  a  change  in  its 
position.  It  should  be  noted  that,  fluctuations  in  other  static  quantities,  such  as  the  density  and 
temperature,  can  also  introduce  perturbations  to  the  equilibrium  conditions,  (i.e.,  the  balance  of 
the  downstream  and  upstream  conditions  across  the  separation  shock).  Note  that  this  is  strictly 
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a  local  phenomenon  and  fluctuations  in  static  quantities  must  travel  to  the  vicinity  of  the  separa¬ 
tion  shock  to  generate  perturbations  in  the  equilibrium  conditions.  In  these  flows,  fluctuations 
are  generally  caused  by  turbulence,  and  can  either  propagate  in  the  flow  or  be  convected  by 
turbulent  structures  to  the  vicinity  of  the  separation  shock.  It  should  also  be  noted  that,  in  sep¬ 
arated  interactions,  a  substantial  part  of  the  flow  downstream  of  the  shock  can  be  subsonic  and 
turbulent  fluctuations  can  propagate  or  be  convected  in  the  upstream  direction.  These  turbulent 
fluctuations  generally  occur  at  high  frequencies  and,  therefore,  can  only  induce  high  frequency 
fluctuations  in  the  ratio  of  static  quantities  and  thus  the  separation  shock  velocity.  This  part  of 
the  mechanism  is  identical  in  all  interactions. 

The  second  part  of  the  mechanism  is  related  to  the  dynamics  of  the  separation  bubble. 
In  this  part  of  the  model,  the  separation  bubble  is  considered  as  a  “flexible”  and  “permeable” 
body  of  fluid  which  behaves  as  an  extension  of  the  shock  generating  geometry.  The  separation 
shock  forms  at  the  upstream  end  of  the  bubble.  The  expansion  of  the  separation  bubble,  for 
example,  displaces  the  separation  shock  to  an  upstream  position,  whereas  the  contraction  of  the 
bubble  results  in  a  displacement  in  the  downstream  direction.  In  addition  to  this,  the  trembling 
of  the  separation  bubble  can  also  cause  changes  in  separation  shock  position.  The  characteris¬ 
tic  frequencies  of  this  unsteadiness  of  the  separated  flow  are  different  in  different  interactions. 
Among  other  things,  they  depend  on  the  size  and  the  internal  structure  and  dynanucs  of  the 
separated  flow.  In  general  and  in  a  given  interaction,  the  larger  the  separated  flow  the  lower 
the  characteristic  frequencies.  However,  the  unsteadiness  of  a  certain  size  separated  flow  from 
the  UCR  interaction,  for  example,  is  characterized  by  much  lower  frequencies  than  that  of  a 
comparable  size  separated  flow  from  the  SCR-25  interaction.  This  is  because  the  internal 
structure  and  dynamics  of  the  latter  are  quite  different.  In  the  SCR-25  interaction,  the  separated 
flow  is  actually  a  spiraling  vortex  rather  than  the  nominally  two-dimensional  separation  bubble 
in  the  UCR  interaction.  The  unsteadiness  of  the  spiraling  vortex  is  mostly  in  the  forin  of  a 
trembling  type  motion  occurring  at  high  frequencies,  whereas  that  of  the  latter  is  usually  in  the 
form  of  an  expansion  and  contraction  type  motion  occurring  at  much  lower  frequencies.  The 
unsteadiness  of  the  separated  flow  may  be  caused  by  such  physical  phenomena  as  the  entrain¬ 
ment  of  turbulent  structures,  which  would  expand  the  bubble,  or  it  could  simply  be  a  dynamic 
response  of  this  “flexible”  body  of  fluid  to  incoming  boundary  conditions.  In  the  latter,  a  fuller 
incoming  velocity  profile  would  cause  the  separation  to  occur  further  downstream  leading  to  a 
smaller  separation  bubble  or  a  displaced  bubble  of  the  same  size.  Thus  variations  in  the  veloc¬ 
ity  profile  of  the  incoming  turbulent  boundary  layer  can  result  in  a  different  size  separation 
bubble  or  in  a  “sliding”  motion  for  the  same  size  bubble.  It  should  be  noted  that,  although  the 
physical  mechanism  describing  the  large-scale,  global  motion  of  the  separation  shock  is  the 
same  for  all  the  interactions  (i.e.,  expansion  and  contraction  of  the  bubble  and/or  its  trembling), 
this  mechanism  is  characterized  by  different  frequencies  in  different  interactions. 


4.8.  Source  of  low-frequency  component 

From  the  work  described  above,  it  is  clear  that  the  high-frequency,  jittery  motion  of  the 
separation  shock  is  the  result  of  the  passage  through  the  wave  of  individual  large-scale  turbu¬ 
lent  structures.  Thus,  it  is  reasonable  to  conclude- that  this  component  of  the  unsteadiness  is  an 
inherent  feature  of  all  turbulent  flows.  The  primary  outstanding  question  concerns  the  cause  of 
the  low-frequency  expansion/contraction  of  the  separated  flow  which  is  characterized  by  the 
large-scale,  long-duration  excursions  of  the  separation  shock  wave. 

Preliminary  experimental  work  to  address  this  question  has  revealed  two  very  interest¬ 
ing,  complementary  results.  First,  there  is  a  distinct  correlation  between  large-scale  exppsion 
or  contraction  of  the  separated  flow  and  long  duration  (i.e.,  low-frequency)  falls  or  rises  m 
pitot  pressure  in  the  incoming  turbulent  boundary  layer.  In  these  experiments,  fluctuating  pitot 
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pressures  were  measured  9.28  inches  (13.35 7 ,  where  5/  =  mass  flux  thickness)  upstream  of 
a  compression  ramp  comer,  while,  at  the  same  time,  fluctuating  wall  pressures  were  measured 
within  the  intermittent  region.  Specific  motions  of  the  shock  foot  were  then  used  as  “triggers” 
for  extraction  of  ensembles  from  the  fluctuating  pitot  signal.  Ensemble-average  pitot  pressure 
signals  are  presented  in  Fig.  14a  for  an  upstream  shock  sweep  and  in  Fig.  14b  for  a  down¬ 
stream  shock  sweep.  Within  the  boundary  layer  (F  <  0.95; ),  there  is  a  well  resolved  decrease 
in  pitot  pressure  while  the  shock  is  undergoing  an  upstream  motion,  and  an  increase  in  pitot 
pressure  while  the  shock  is  undergoing  a  downsteam  motion.  The  time  scale  of  these  gradual 
rises  and  falls  in  pitot  pressure  is  around  1  ms,  almost  two  orders  of  magnitude  higher  than  the 
large  eddy  time  scale  S 7  /C/«,.  Thus,  these  rises  and  falls  in  pitot  pressure  are  not  generated  by 
individual  turbulent  structures,  but  are  caused  by  low-frequency  “global”  motion  of  the  bound¬ 
ary  layer.  Second,  results  from  the  same  experiment  show  that  the  ensemble-averaged  pitot 
pressure  at  a  fixed  location  in  the  incoming  undisturbed  boundary  layer  correlates  with  separa¬ 
tion  shock  wave  position.  Three  positions  were  considered;  (1)  shock  foot  at  upstream  edge 
of  intermittent  region  (i.e.,  corresponding  to  large  separation  bubble),  (2)  shock  foot  at 
downstream  edge  of  intermittent  region  (i.e.,  small  separation  bubble),  and  (3)  shock  foot 
centered  within  intermittent  region.  Results  showed  that  the  ensemble-average  pitot  pressure 
was  below  its  overall  mean  value  for  case  1  (large  bubble),  above  it  for  case  2  (small  bubble), 

and  equal  to  it  for  case  3.  ,  .  ,  j 

These  two  sets  of  results  suggest  a  model  in  which  the  boundary  layer  thickens  and 
thins”  as  it  flows  downstream  (that  is  why,  at  a  fixed  point  above  the  wall,  the  pitot  pressure 
rises  or  falls).  In  such  a  simple  model,  a  higher  pitot  pressure  at  a  fixed  height  would  imply  a 
fuller  velocity  profile  than  the  mean  profile,  one  more  resistant  to  separation,  which  would 
generate  a  smaller  separation  bubble  than  would  the  mean  profile.  Conversely,  a  lower  pitot 
pressure  at  a  fixed  height  would  imply  a  slightly  retarded  profile  which  would  generate  a  larger 
separation  bubble  than  would  the  mean  profile.  Both  observations  are  consistent  with  the 
experimental  observations,  as  is  the  fact  that,  when  the  shock  foot  is  at  its  mean  position  (case 
3  above),  the  ensemble-averaged  pitot  pressure  takes  the  overall  mean  value.  However,  at  this 
stage  it  must  be  emphasized  that  the  “thickening-thinning”  boundary-layer  model  is  an 
inference. 

In  speculating  about  possible  causes  of  a  thickening-thinning  boundary  layer,  the 
question  of  Taylor-Gortler  vortices  arose.  When  a  boundary  layer  develops  on  a  concave  wall, 
Gortler  vortices  can  occur  within  the  boundary  layer  (Fig.  15).  A  centrifugal  force  imposed  by 
the  curvature  of  the  wall  forces  fluid  particles  to  be  displaced  towards  the  wall.  On  the  other 
hand  the  pressure  force  opposes  this  centrifugal  force  to  restore  the  particles  to  their  original 
position.  When  the  centrifugal  force  is  larger  than  the  pressure  force  such  that  the  viscosity 
cannot  damp  the  motion,  the  flowfield  is  considered  to  be  unstable  and  small  disturbances  are 
likely  to  be  amplified.  This  is  the  underlying  mechanism  of  the  counter-rotating  streamwise 
vortices,  usually  called  “Gortler  vortices.”  Although  Gortler  instability,  which  is  inviscid  in 
nature,  has  undergone  many  investigations,  these  investigations  are  mostly  limited  to  laminar 
and  incompressible  flowfields  with  emphasis  on  their  role  in  the  transition  process.  The  exis¬ 
tence  of  these  vortices  in  compressible,  turbulent  flowfields,  and  how  far  downstream  of  a 
region  of  concave  curvature  such  vortices  persist  are  less  clear.  For  a  detailed  overview  of  cur¬ 
rent  knowledge  of  the  subject,  readers  are  referred  to  the  excellent  recent  reviews  by  Floryan 
and  Saric  (1994). 

Earlier  work  (Hoffman  et  al.  1985)  shows  that  vortices  developing  on  a  concave  wall 
wander  about  preferred  lateral  positions  and  these  positions  are  likely  to  be  determined  by 
nonuniformities  in  the  oncoming  flow.  Although  the  spanwise  wavelength  is  usually  about 
twice  the  boundary  layer  thickness,  the  spacing  and  amplitude  are  functions  of  the  initial 
disturbances.  In  his  review,  Floryan  suggests  that  the  nonuniformities  in  the  boundary  layer 
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lead  to  spanwise  variations  in  time-averaged  measurements  which  signal  the  presence  of  vor¬ 
tices.  Although  this  opinion  is  widely  accepted,  the  streamwise  evolution  of  the  vortices  is  a 
point  of  controversy.  While  some  researchers  (So  et  al.  1975;  Prabhu  et  al.  1981)  observe 
that  these  vortices  decay,  others  (Hoffman  et  al  1985;  Smits  et  al  1981)  do  not  confirm  this 
observation.  Floryan  states  that  a  study  by  Smits  et  al  suggests  that  the  vortices  generated  in 
the  curved  section  of  a  wind  tunnel  persist  almost  indefinitely  downstream  in  the  straight  sec¬ 
tion.  They  observed  large  spanwise  variations  in  the  time-averaged  skin  friction,  streamwise 
and  vertical  velocity  components  and  Reynolds  stresses  as  well  as  the  spanwise  wandering  of 
these  quantities.  In  another  study  at  Mach  5  with  two  different  nozzle  geometries,  Beckwith 
and  Holley  (1981)  investigated  the  onset  of  transition  under  extremely  low  incident  disturbance 
levels.  Their  oil  flow  visualization  results  show  clear  streak  patterns  which  indicate  the  exis¬ 
tence  of  the  Gortler  vortices  up  to  a  certain  unit  Reynolds  number  which  depends  on  the 
geometry  of  the  nozzle.  Above  this  Reynolds  number  they  saw  fine  streak  patterns  which 
could  be  evidence  of  the  vortices  in  turbulent  flow  and  commented  that  either  Gortler  vortices 
persist  far  downstream  in  the  turbulent  flow  or  that  the  boundary  layer  there  has  some  quasi- 
turbulent  state. 

A  well-known  effect  of  Gortler  vortices  is  alteration  of  the  velocity  field  (Fig.  16).  The 
low  speed  fluid  close  to  the  wall  will  be  carried  up  into  a  region  of  high  velocity  and  conse¬ 
quently  will  cause  a  local  decrease  in  the  velocity  of  the  flow  (upwash).  If  the  vortex  pair 
rotates  in  the  opposite  direction,  then  the  high  speed  fluid  will  be  carried  down  to  a  low-speed 
region  and  will  cause  a  local  increase  in  the  flow  velocity  (downwash).  A  common  way  to 
identify  the  existence  of  Gortler  vortices  is  through  surface  flow  patterns.  Usually,  a  dark  sur¬ 
face  is  covered  with  a  thin  layer  of  light-colored  tracer  material.  The  locations  corresponding  to 
a  downward  motion  induced  by  the  vortex  pairs  are  determined  by  the  black  streaks  which  in¬ 
dicate  an  increase  in  shear  stress.  Conversely,  light-colored  streaks  are  indications  of  an  up¬ 
ward  motion  and  a  decrease  in  shear  stress.  The  spanwise  distance  between  two  neighboring 
streaks  of  the  same  color  corresponds  to  one  pair  of  counter-rotating  vortices  and  defines  the 
wavelength  of  the  vortex  system. 

To  understand  whether  such  a  vortex  system  might  exist  in  the  current  facility,  an 
experiment  was  conducted  in  which  two  transducers  with  fixed  spacing  were  located  at  differ¬ 
ent  spanwise  locations.  In  the  absence  of  any  vortex  system  the  cross  correlation  coefficients 
obtained  from  different  pairs  of  transducers  should  be  the  same.  If  there  is  a  vortex  system 
then  it  is  possible  that  differences  would  be  observed  at  different  spanwise  locations  on  the 
wall.  Fig.  17a  shows  the  maximum  cross  correlation  coefficient  as  a  function  of  the  spanwise 
location.  Note  that  Z=3  inches  is  the  centerline  of  the  tunnel.  There  is  clearly  a  systematic 
spanwise  variation  in  the  maximum  cross  correlation  coefficient.  There  is  strong  evidence  of 
periodicity  with  a  wavelength  of  approximately  twice  the  boundary  layer  thickness.  This 
result  suggests  that  Gortler  vortices  are  present  in  the  facility.  A  hypothetical  vortex  pair  sys¬ 
tem  based  on  the  variations  in  the  maximum  cross  correlation  coefficient  is  also  shown  in  Fig. 
17b.  Note  that  the  local  maximum  value  of  R[0,AZ,0lT]max  corresponds  to  the  region  in 
which  the  low  speed  fluid  is  carried  up  into  a  high  speed  flow  region.  This  result  may  implic¬ 
itly  be  related  to  the  variations  in  the  wall  shear  stress  ;  when  the  shear  stress  is  smaller,  this 
will  cause  a  decrease  in  the  friction  velocity,  I/^.  If  inner  scaling  is  used  to  plot  the  maximum 
spanwise  cross  correlation  coefficient,  it  is  evident  that  a  decrease  in  will  effectively  cause 
a  decrease  in  the  dimensionless  parameter  AZt/^  /  for  a  fixed  spanwise  separation.  This, 
in  turn,  will  cause  an  effective  increase  in  the  maximum  cross  correlation  coefficient. 

Examination  of  the  spanwise  behavior  of  the  rms  of  the  fluctuations,  and  the  skewness 
and  flatness  coefficients  (Fig.  18a,  18b  and  18c,  respectively),  does  not  reveal  any  obvious 
periodicity.  As  noted  earlier  there  is  quite  a  large  uncertainty  in  the  rms  measurement  which 
may  mask  small  spanwise  variations. 
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The  spanwise  variations  in  the  displacement  thickness,  5  ,  and  the  skin  friction  coeffi¬ 
cient,  Cf,  obtained  by  Barter  and  Dolling  (1993),  are  given  in  Fig.  19.  Although  it  is  difficult 
to  reach  a  firm  conclusion,  it  can  be  seen  that  there  is  a  local  maximum  in  S  and  local  mini¬ 
mum  in  Cf  very  close  to  the  center  of  the  tunnel.  The  local  minimum  for  5  is  approximately 
0.7  inches  to  the  left  and  right  of  the  local  maximum  close  to  the  center.  This  is  followed  by 
local  maxima  which  yields  a  periodicity  of  approximately  2  5q  .  Similar  remarks  can  be  made 
about  the  skin  friction  coefficient. 

From  the  behavior  of  the  maximum  cross  correlation,  displacement  thickness  and  skin 
friction  coefficient,  it  might  be  inferred  that  the  Gortler  vortices  exist  in  the  test  section.  How 
these  vortices  affect  the  dynamics  of  shock-induced  separation,  if  at  all,  is  an  interesting  ques¬ 
tion.  For  close  spacings,  the  spanwise  cross  correlation  of  wall  pressure  signals  has  the 
expected  shape  (Fig.  20),  the  time  delay  at  maximum  correlation  is  zero,  and  the  decay  of  the 
maximum  correlation  coefficient  is  rapid  at  least  for  AZ  /  5o^l.  All  of  this  is  consistent  with 
pressure  signals  generated  by  large  scale  turbulent  structures.  Beyond  AZ  I  5q  of  about  1, 
the  decay  is  more  gradual  and  at  certain  spacings  unusual  cross  correlations  are  observed  (i.e.. 
Fig.  21).  Non-zero  time  delays  such  as  those  in  Fig.  21  must  be  caused  by  structures,  or  other 
pressure  generating  mechanisms,  oblique  to  the  flow  direction.  No  proof  can  be  offered  at  this 
stage  but  the  question  arises  as  whether  the  gradual  decay  of  the  cross  correlation  coefficient 
and  results  of  Fig.  21  could  be  related  to  the  spanwise  wandering  of  Gortler  vortices  about 
their  preferred  positions,  which  in  turn  also  affect  the  shock  foot  motion. 


4.9.  Results  from  computational  analysis 

To  examine  the  validity  of  the  first  part  of  the  physical  model  developed  from  experi¬ 
mental  results,  a  computational  analysis  has  also  been  performed.  In  this  study,  time-accurate 
Euler  equations  governing  the  inviscid  flow  in  a  quasi  one-dimensional  channel  have  been 
solved  using  the  Van  Leer  flux- vector  splitting  scheme  (Van  Leer  et  al.  1987).  The  inlet  and 
outlet  boundary  conditions  were  so  prescribed  to  cause  the  formation  of  a  normal  shock  in  the 
middle  of  the  channel  for  the  steady-state  solution.  The  channel  geometry,  and  the  Mach  num¬ 
ber  and  static  pressure  distributions  for  the  steady-state  solution  are  plotted  in  Figs.  22a,  22b, 
and  22c,  respectively.  Since  the  convergence  criterion  was  based  on  a  global  average-epor 
calculation,  the  shock  is  smeared  over  a  couple  of  cells.  For  the  purposes  of  this  study,  this  is 
not  important.  Once  the  steady-state  solution  converged,  the  static  pressures  at  the  inlet  and 
outlet  boundary  conditions  were  varied  sinusoidally  using  the  summation  of  four  different 
sinusoids.  This  permits  different  perturbations  by  changing  the  amplitude  and  frequency  of 
each  of  the  four  sinusoids.  Typical  pressure  disturbance  histories  introduced  at  the  inlet  and 
outlet  are  plotted  in  Figs.  23a  and  23b,  respectively.  The  response  of  the  normal  shock  to 
these  perturbations  was  monitored  by  calculating  its  position  (i.e.,  the  cell  in  which  it  occurs) 
in  the  channel  at  each  time  step.  Since  the  size  of  the  cells  dictates  the  spatial  resolution  by 
which  the  shock  position  can  be  determined,  a  simple  linear  interpolation  has  been  used  to  yield 
a  piecewise  smooth  normal  shock  position  history,  Xc(t).  The  normal  shock  velocity  history, 
Vc(t)  can  be  calculated  by  simply  taking  the  temporal  derivative  of  its  position  history,  Xc(t)- 
The  normal  shock  position  and  velocity  histories,  corresponding  to  the  disturbances  shown  in 
Figs.  23a  and  23b,  are  plotted  in  Figs.  23c  and  23d,  respectively.  The  static  pressure  ratio 
across  the  normal  shock,  Rc(t)  is  also  calculated  at  each  time  step.  If  the  normal  shock  is  in  cell 
i,  for  example,  then  the  upstream,  pu(t),  and  downstream,  pd(t),  static  pressures  are  calculated 
from  the  conservative  variables  in  cells  i-2,  and  i+2,  respectively.  The  data  from  cells  imme¬ 
diately  upstream  {i-1)  and  downstream  (i+i)  of  the  shock  were  not  used  for  this  calculation. 
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since  the  shock  is  smeared  over  a  couple  of  cells  as  mentioned  earlier.  The  static  pressure  ratio 
history  generated  by  the  inlet  and  outlet  disturbances  is  plotted  in  Fig.  23e. 

Once  sufficient  amounts  of  data  for  the  normal  shock  velocity,  Vc(t)  and  static  pressure 
ratio,  Rc(t)  fluctuations  were  collected,  the  process  was  stopped.  These  data  were  then  corre¬ 
lated  using  the  variable  interval  cross-correlation  method  described  earlier.  The  normalized 
results  of  this  correlation  are  shown  in  Fig.  24.  This  correlation  has  a  very  strong  negative 
peak  which  occurs  at  a  negative  time  delay.  This  shows  very  clearly  that  fluctuations  in  the 
static  pressure  ratio  across  the  normal  shock  cause  opposite  fluctuations  in  its  velocity  (also  see 
Figs.  23d  and  23e).  Furthermore,  the  very  strong  correlation  coefficient  suggests  that  this  is 
the  only  physical  mechanism  responsible  for  the  normal  shock  oscillations  in  the  channel. 

Even  though  these  results  are  from  a  quasi  one-dimensional  analysis,  they  nevertheless 
demonstrate  the  unsteadiness  caused  by  fluctuations  in  the  ratio  of  static  quantities  as  dictated 
by  the  governing  equations.  This  analysis,  therefore,  provides  computational  support  for  the 
first  part  of  the  physical  mechanism  developed  from  experimental  results,  as  can  be  seen  by 
comparing  the  results  presented  in  Fig.  5  with  those  in  Fig.  24. 


5.  Conclusions 

The  purpose  of  this  work  was  (i)  to  examine  separation  shock  unsteadiness  in  different 
interactions  and  determine  whether  a  universal  model  describing  the  unsteadiness  could  be 
developed,  and  (ii)  to  determine  whether  or  not  the  observed  unsteadiness  is  a  feature  of  turbu¬ 
lent  flows  in  general,  or  is  specific  to  the  wind  tunnel  environment.  To  this  end,  wall  pressure 
and  pitot  pressure  fluctuation  measurements  were  made  in  interactions  generated  by  unswept 
and  25  deg  swept  compression  ramp  models,  and  by  8  deg  and  30  deg  swept  blunt-fin  models 
in  a  high  Reynolds  number,  Mach  5  turbulent  boundary  layer.  Using  a  new  conditional  sam¬ 
pling  algorithm,  the  separation  shock  unsteadiness  was  characterized  in  terms  of  its  position, 
Xs(t)  and  velocity,  Vs(t)  histories.  These  fluctuations  were  then  correlated  with  the  condition¬ 
ally  extracted  static  pressure  ratio,  Rs(t)  across  the  separation  shock,  and  with  wall  pressure 
measurements  made  upstream  and  downstream  of  the  intermittent  region  in  each  of  these  inter¬ 
actions.  A  comprehensive  physical  model  describing  the  unsteadiness  has  been  developed.  In 
this  model,  the  separation  shock  unsteadiness  can  be  caused  by  two  different  physical 
phenomena  which  comprise  the  two  parts  of  the  mechanism. 

In  the  first  part,  the  small-scale  or  jittery  motion  of  the  separation  shock  is  caused  by 
fluctuations  in  the  ratio  of  static  quantities  across  the  shock  foot.  These  perturbations  induce 
fluctuations  in  the  separation  shock  velocity  and  thus  changes  in  separation  shock  position  are 
an  integral  result  of  the  velocity  fluctuations.  An  increase  in  the  ratio  of  static  pressure,  for  ex¬ 
ample,  causes  an  upstream  excursion  in  the  separation  shock  position  and  vice  versa.  The  ratio 
of  static  quantities  can  be  changed  by  independent  fluctuations  in  upstream  and  downstream 
conditions  or  both.  These  fluctuations  are  attributed  to  turbulence  and  can  either  propagate  or 
be  convected  to  the  vicinity  of  the  separation  shock. 

In  the  second  part,  the  large-scale  or  global  motion  of  the  separation  shock  is  caused  by 
the  expansion  and  contraction  or  trembling  motion  of  the  separation  bubble.  Here,  the  separa¬ 
tion  bubble  is  considered  as  a  “flexible”  and  “permeable”  body  of  fluid  which  behaves  as  an 
extension  of  the  shock  generating  geometry.  The  separation  shock  forms  at  the  upstream  end 
of  the  bubble.  Thus,  an  expansion  of  the  separation  bubble,  for  example,  displaces  the  separa¬ 
tion  shock  to  an  upstream  position,  whereas  a  contraction  of  the  bubble  results  in  a  displace¬ 
ment  in  the  downstream  direction.  The  unsteadiness  of  the  separation  bubble  may  be  caused 
by  such  physical  phenomena  as  the  entrainment  of  turbulent  structures,  which  would  expand 
the  bubble,  or  it  could  simply  be  a  dynamic  response  of  this  “flexible”  body  of  fluid  to  incom¬ 
ing  boundary  conditions.  Thus  variations  in  the  velocity  profile  of  the  incoming  turbulent 
boundary  layer  can  result  in  a  different  size  separation  bubble  or  in  a  “sliding”  motion  for  the 
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same  size  bubble.  It  should  be  noted  that  some  evidence  has  been  found  of  Taylor-Gortler  vor¬ 
tices  in  the  incoming  flow,  and  it  is  possible  that  these  play  a  role  in  the  low-frequency 
physics.  Confirmation  using  intrusive  instrumentation  alone  is  extremely  difficult,  and,  in  the 
next  phase  of  the  work,  both  conventional  instrumentation  and  planar  laser  imaging  methods 
will  be  used  to  explore  the  incoming  flow  structure  corresponding  to  different  separated  bubble 
sizes  and  dynamics. 

In  addition  to  experimental  results,  computations  in  an  inviscid,  unsteady  flow  have 
also  been  performed.  In  this  study,  the  response  of  a  normal  shock  in  a  quasi  one-dimensional 
channel  was  examined  as  perturbations  in  upstream  and  downstream  boundary  conditions  were 
introduced.  Although  this  does  not  model  the  experimental  conditions,  the  results  clearly 
demonstrate  the  unsteadiness  caused  by  fluctuations  in  the  ratio  of  static  quantities  as  dictated 
by  the  governing  equations.  This  analysis,  therefore,  provides  computational  support  for  the 
first  part  of  the  physical  mechanism  developed  from  experimental  results. 
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Fig.  1  Unswept  and  Swept 


Fig.  3  Normalized  Distributions  of  Mean, 
Standard  Deviation,  and  Intermittency 


Fig.  4  Normalized  Power  Spectral  Densities  for 
Separation  Shock  Position,  Xs(t),  and  Velocity, 
Vs(t)  Histories 


Fig.  5  Normalized  Cross-Correlation  of  Shock 
Velocity,  Vs(t),  and  Pressure  Ratio,  Rs(t) 
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Fig.  6  Normalized  Cross-Correlation  of  Shock 
Velocity,  VsCt),  and  Upstream  Pressure,  P^CO 


Fig.  7  Normalized  Cross-Correlation  of  Upstream 
and  Downstream  Pressure  Measurements 


Fig.  8  Normalized  Cross-Correlation  of  Shock 
Velocity,  Vs(t),  and  Downstream  Pressure, 


Fig.  9  Normalized  Cross-Correlation  of  Shock 
Velocity,  Vs(t),  and  Downstream  Pressure, 
(High-Pass  Filtered,  3  <  f  <  50  kHz) 


Fig.  10  Normalized  Cross-COTrelation  of 
Downstream  Pressure  Measurements 
(High-Pass  FUtered,  3  <  f  <  50  kHz) 


Fig.  1 1  Normalized  Cross-Correlation  of  Shock 
Velocity,  Vs(t),  and  Downstream  Pressure,  P^ft) 
(Low-Pass  Filtered,  0  <  f  <  3  kHz) 
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Fig.  12  Normalized  Distributions  of  Ensemble- 
Averaged  Pressure  Measurements  (Erengil  & 
Dolling  1991a) 


Fig.  13  Normalized  Cross-Correlation  of  Shock 
Position,  Xs(t),  and  Downstream  Pressure,  P^Ct) 
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Figure  14a.  Ensemble-averaged  pitot  pressure  signals  in  the 
incoming  boundary  layer  for  an  upstream  shock  sweep. 
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Figure  14b.  Ensemble-averaged  pitot  pressure  signals  in  the 
incoming  boundary  layer  for  a  downstream  shock  sweep. 


Figure  17.  a)  Variation  of  maximum  cross-correlation 
coefficient  across  tunnel  floor  span;  b)  hypothetical  vortex 
system. 


Figure  18.  Spanwise  variation  of  a)  rms  wall  pressure; 
b)  skewness  coefficient;  c)  flatness  coefficient. 


Figure  19.  Spanwise  variation  in  displacement  thickness  and 
in  skin  friction  coefficient. 


AZ/5o 


Figure  20.  Decay  of  maximum  cross-correlation  coefficient  as 
a  function  of  Az/ Sq- 
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Figure  21.  Example  of  span  wise  cross  correlation  with 
unexpected  features. 


Figure  22.  Channel  geometry,  and  steady-state  solutions  for 
Mach  number  and  pressure. 
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velocity,  Vc(t)  and  pressure  ratio,  Rc(t). 


